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Abstract 

Motivated by the desire to test modified gravity theories exhibiting the Vainshtein mechanism, we solve in various 
physically relevant limits, the retarded Galileon Green's function (for the cubic theory) about a background 
sourced by a massive spherically symmetric static body. The static limit of our result will aid us, in a forthcoming 
paper, in understanding the impact of Galileon fields on the problem of motion in the solar system. In this 
paper, we employ this retarded Green's function to investigate the emission of Galileon radiation generated by 
the motion of matter lying deep within the Vainshtein radius r v of the central object: acoustic waves vibrating 
on its surface, and the motion of compact bodies gravitationally bound to it. If A is the typical wavelength of 
the emitted radiation, and r is the typical distance of the source from the central mass, with r -C r v , then, 
compared to its non- interacting massless scalar counterpart, we find that the Galileon radiation rate is suppressed 
by the ratio (7\,/A)~ 3 / 2 at the monopole and dipole orders at high frequencies r„/A » 1. However, at high enough 
multipole order, the radiation rate is enhanced by powers of r v /r a . At low frequencies r v /X <C 1, and when the 
motion is non-relativistic, Galileon waves yield a comparable rate for the monopole and dipole terms, and are 
amplified by powers of the ratio r v /r for the higher multipoles. 



I. INTRODUCTION AND MOTIVATION 



The study of the problem of motion in General Relativity (GR) has been a central theme in testing 
its validity. For instance, the post-Newtonian program, computing GR corrections to the Newtonian 
gravitational potential between massive bodies, is crucial to understanding gravity in our solar system. 
In the past decade or more, the post-Newtonian analysis of weak field gravity has also been developed 
to very high order in perturbation theory because of the need to model gravitational waves (GWs) 
from inspiraling compact binaries, which is expected to be a major source for detectors like Advanced 
LIGO. At the same time, the discovery of cosmic acceleration and its associated cosmological constant 
problem, has prompted many attempts to modify how gravity operates at large (astrophysical) length 
scales. One such example is the family of scalar field theories known as Galileons [I], building on 
interesting properties of a limit of the Dvali-Gabadadze-Porrati model These scalar fields couple 
to the stress-energy of ordinary matter, and therefore alter large scale dynamics, but due to their self 
interactions they exhibit the Vainshtein screening effect (3j 0], such that masses close to the central 
source of gravity (lying within a so-called Vainshtein radius r v ) do not feel their presence, and thus 
allowing Galileons to evade solar system tests of GR to date. Furthermore, models of this type have 
become even more interesting since it has been discovered that they may be extended in multiple ways 
to yield other new field theories [5H24| with related attractive properties, and that they arise as a limit 
of ghost-free massive gravity [25H27j 

In this paper, we wish to lay the groundwork for understanding analytically, and in some detail, the 
impact of such Vainshtein screened scalar fields on the problem of motion. To achieve concrete results 
we will consider the cubic Galileon theory about flat spacetime and couple the Galileon field to the trace 
of the stress-energy of matter. We will place a large mass M at the origin of our coordinate system. The 
central goal of this paper is to solve the retarded Green's function of the linearized equations of motion 
of the Galileon fluctuating around the background field sourced by M. In a paper in preparation [28], we 
use the static limit of our Green's function here to investigate the conservative portion of the dynamics; 
to compute the effective potential between well separated test masses orbiting around M. 

In the current paper we address the dissipative aspect of the dynamics: does motion of matter lying 
well within the Vainshtein radius of M produce Galileon radiation that can carry energy-momentum 
away to infinity? This question arises as an issue of principle because the Galileon is a massless scalar 
field, and one would expect the motion of sources of massless fields to create radiation. Yet, it is not 
clear how much radiation would actually be produced, because it could perhaps be suppressed by the 
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Vainshtein mechanism. 

Beyond issues of principle, we are also motivated by the possibility that one can constrain modified 
gravity theories by demanding that the power loss from binary pulsar systems, such as the Hulse- Taylor 
binary PSR B1913+16, had better not deviate too far from the predictions of GR, since observations 
have confirmed the latter to high precision 1 . Furthermore, we look forward to the prospect that, as 
already alluded to, within the next decade or so GW detectors may be able to directly listen in on 
the spacetime ripples generated by such compact binary systems. Once this is possible, we may hope 
to use these GW signals to search for or put further constraints on the existence of the Vainshtein 
mechanism. This requires that we develop a quantitative prediction of the Galileon waves themselves, 
beyond just an estimate of the power loss through scalar emission. Even though the two comparable 
point mass Galileon problem possibly requires different techniques to solve due to the importance of the 
nonlinearities of the field equations, the results of our current paper, which assume the existence of a 
very large central mass, may perhaps be seen as an approximation to the situation where the inspiraling 
binary consists of rather unequal masses, say Mi <C M2. Yet another possible source of GWs is the 
oscillations of neutron stars themselves; as such, we will also consider a toy problem of surface waves 
on a spherical body stimulating Galileon waves. 

The outline of the paper is as follows. In section (fill), we set up our problem in a quantitative 



manner. In section (III A) we summarize the results for the Galileon retarded Green's function about 



the background field of the central mass M; and following that in section (IIIB) we step through its 



derivation. In section (IIIC) we describe how the linearized equations of motion of the Galileon field 



propagating about the background sourced by M is equivalent to a minimally coupled massless scalar in 



some curved spacetime. Then in section (IV) we move on to use the retarded Galileon Green's function 
to study two examples of radiative processes: Galileon radiation produced by surface waves on an 
otherwise spherical central body, and that generated by n point masses gravitationally bound to M. In 
appendix (|A|), we work out both the curved and Minkowski spacetime minimally coupled massless scalar 
analog of the Lienard-Wiechert potentials in electromagnetism; this spacetime calculation complements 



the frequency space one in section ( IV B ) . 

Notation A few words on notation. The speed of light is set to unity. The Galileon field is 
II. The background spacetime is (3+l)-dimensional Minkowski, with metric in Cartesian coordinates 



For recent work on this topic, see 
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given by 

ds 2 = r]^dx^dx u , 7]^ = diag[l, -1,-1,-1] (1) 

so that, 

d 2 u = v^d^u, (an) 2 = r^^n^n. (2) 

In spherical coordinates (t, r, 9, cf>), with £ [0, tt] and G [0, 2ir), the metric reads instead 

ds 2 = dt 2 - dr 2 - r 2 n AB dx A dx B , (3) 
^AB = diag[l,sin 2 0], x A = (e,(f ) ), (4) 

where df2 = d6d(j>\f& (with \/f2 = V det S7ab = sin #) is the infinitesimal solid angle. 

A hat on a variable representing spatial location, e.g. x, denotes the unit vector, x = x/\x\. In 
particular, x only depends on the spherical coordinate angles x = x[8,(/)]. Finally, the Planck mass is 
defined in terms of Newton's gravitational constant Gn as M p \ = 1/y / 327tGn. 



II. SETUP 



We would like to understand Galileon radiative processes taking place in the background Galileon 
field II[r] generated by a massive central body of mass M, which we would take to be static (time 
independent) and spherically symmetric. To model this mass M we shall treat it as a point particle at 
rest, located at the spatial origin of the coordinate system. The Galileon radiation we are investigating 
is generated by matter (described by stress-energy tensor 5T^ U ) lying well within the Vainshtein radius 
r v of the central mass. (The Vainshtein radius r v , as we will see very shortly, is the radius below which 
the Galileon field II generated by the mass M is increasingly governed by non- linear self interactions; 
well outside Vainshtein, the theory is linear and there M generates a 1/r Coulomb potential.) This 
matter distribution bT^ v is meant to be viewed as a perturbation relative to the mass M, but can 
otherwise be arbitrary. For instance 5T^ V may describe slight deviation of the mass M from an exact 
spherical configuration (its multipole moments), and/or a deviation from time independence; it could 
also be the energy-momentum of n light compact bodies, with masses {m a <C M\a = 1, . . . , n}, orbiting 
around M. 

We will take for the Galileon theory the simplest cubic II Lagrangian; in actuality, several other 
terms are allowed in 4 spacetime dimensions. The Galileon, being a scalar, couples to the trace of the 
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stress-energy tensor of matter. The total action for our setup is therefore 

Sn + S M + 5S (5) 

where, in Cartesian coordinates, 

s n = J a'x Q(an) 2 + ia 2 n(9n) 2 ) , a > o, (6) 
Sm^^Ju'WA, (7) 

SS = I d 4 x-^<5T, 5T = 5T» (8) 
J M p i " v ' 

Because we are interested in radiation that can propagate to infinity, we note that far away from the 
source M, contributions to the Galileon stress-energy tensor due to the cubic-in-Il terms in Sit fall 
away more rapidly than their quadratic counterparts (provided, of course, that II and its gradients falls 
off with increasing r). Expressed in Cartesian coordinates, the asymptotic Galileon energy- momentum 
tensor is thus that of the non-interacting massless scalar in flat spacetime, 

-> °°] = dfJIdvIL - ^(dlif . (9) 

The background Galileon field sourced by the mass M is the exact solution to 5{Sj\+Sm) / = 0. Since 
we have a static spherically symmetric source, II must only depend on r, and so the Euler-Lagrange 
equation from varying Sn + Sm becomes 

^ 2 {^-^(a r n) 2 }) =- r ^fM do) 



Integrating both sides of (10) over a sphere of radius r centered at then yields 



_ 4 _ 9 M 

d r u - — (d r uy = - ^ 9 . (ii) 

A 3 r V 1 4vrM pl r 2 v ' 



(Even though we have modeled the central body as a point mass, ( 11 ) is valid outside any isolated static 



spherically symmetric matter distribution as long as we replace M with the integral 4-7T r^[r']r /2 dr', 



where r^ v is the stress-energy tensor of the body.) Eq. (11) is a quadratic equation in II [r]. Defining 
the Vainshtein radius 

1 /4 M \ 1/3 

r » S AU^) ' (12) 



we obtain 

II \r\ r ( I /r v \ 3 \ 

— (13) 



A 3 
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There is a second solution for II [r] in which the negative sign in front of the square root is replaced with 
a plus sign. This solution is proportional to r and hence blows up as r — > oo. Since the stress-tensor 
depends on IT [r] - see eq. ^ and note that including the contributions to from the cubic-in-II 
terms in Sri would only exacerbate the problem - we may discard this second solution on the grounds 
that the energy- momentum of II measured by an asymptotic observer cannot be infinite. 



It is possible to integrate eq. (13) exactly in terms of hypergeometric functions, 

" ( 1 - ! ( 

r 



II[r] 



r 

16 



3 

r v \ 2 



2-^1 



1 1 7 

2' 6' 6' 



r 



r . 



r[-!]r[J] 



where we have chosen the asymptotic boundary condition to be II [r — > 00] = 0. 



For later use, note the following limits of equations (|13|) and (|14|). When r <^ r v , 

2l \ 



H[r] 



M 



n-ir 7 



6 J 



+ 



+ 



and 



n'[r] 



n"[r] 



A 3 ? 




l + O 




+ 2 + 



3 

f \ 2 



On the other hand, when r 3> r v , 



U[r] 



M 



4vrM p ir 



l + O 



and 



n [r] 



n [r] 



A 3 r 
16" 



\ 3 

r 



A 3 



+ 



(14) 



(15) 



(16) 
(17) 



(18) 



(19) 
(20) 



Notice that the first term of the small radius limit II [r] ~ 1 j^fr in ( 16 ) can be obtained by dropping 
the linear term <9 r II[r] in eq. (11), whereas the first term in (19) of the large radius limit Il'fr] ~ 1/r 2 , 
which is the usual Coulomb force law, comes about from dropping the non-linear (<9 r IT) 2 piece in eq. 



(11). This is the quantitative statement that, close to the matter source (r <C r v ), the dynamics of 
Galileons are primarily governed by their nonlinear self interactions. The theory is linear when the 
observer is well outside the Vainshtein radius. 
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More physically, to see the Vainshtein effect at work, imagine a test point mass m <C M, at spatial 
location Z, orbiting the central body. Its action takes the same form as Sm in eq. Q, except we 
evaluate the Galileon field about the background II generated by M, 



m 



M, 



pi 



dm 



Z{t] 



(21) 



(Strictly speaking, Lorentz invariance says dt needs to be replaced with dty 1 — Z[t] 2 ; however, we are 
working in the non-relativistic regime, where the square root is very close to unity.) By employing the 



definition of the Vainshtein radius in eq. (12), followed by integrating the first term on the right-hand- 



sides of equations (16) and (19) without worrying too much about the overall numerical factors, 

G N Mm / \Z\t\ 



\z[t]\ 

G N Mm 



\Z[t]\ <Cr„ 

\Z[t}\ >r„ 



(22) 
(23) 



Thus, the Galileon potential experienced by a test point mass orbiting close to the central mass M is 
the Newtonian gravitational potential G^M/\Z\ multiplied by a suppression factor of {\Z\/r v fl 2 < 1. 
Only when the test mass travels well outside Vainshtein does the suppression factor drop out and the 
Galileon potential become comparable in strength to that of regular gravity. 2 
With the exact background solution Il'[r] in hand, we now substitute 

n[a?] = TT[r] + ip[x] (24) 

in ([5]), and drop terms cubic-in-</?. The resulting linearized equation of motion of ip about the background 
IT reads 

5T[x] 



Wx<p[x] 



M pl ' 



where the differential operator W x is 



1 



W x <p[x] = ( ei9 f 2 - e 2 8 2 r - -e 3 8 r - -^e 3 L 2 ) <p[x] 



(25) 



(26) 



and L 2 is the angular part of the Laplacian in Euclidean 3-space, usually called the negative of the 
"angular momentum squared" operator, given by 



Pip = -^=d A (Vnn^dBv) . 



(27) 



The reader concerned about the stability of the Galileon model is referred to [6]. 
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Here 



ei[r] = 1 

e 2 [r] = 1 
e 3 [r] = 1 



8IT 

srf 

A 3 r 



1 / 3 / (2r3 + rg)" 2 




1 /( lr' + r;^ 



(28) 

(29) 
(30) 



A 3 4 \/ r 3 (r 3 + r 3 ) ' 

If we assume 5T does not implicitly depend on IT, then tp is entirely sourced by 5T, the (for now 
arbitrary) matter perturbations. The key to solving cp in terms of 5T is the retarded Green's function 
defined by the equation, 

1 



W x G[x,x'} = W x >G[x,xf] = —S[t - t']6[r - r'}6[4> - </>']6[cos9 - cos6'], 



rr 



(31) 



where the Ss are the Dirac delta functions. The solution to the linearized Galileon equation ( 25 ) about 
the static spherically symmetric background IT is now (in Cartesian coordinates) 

h 8T[x>] 



tp[x] = J d 4 x'G[x, x'} - 



pi 



(32) 



We write x to represent a collective label for (t, r, 9, <f>) and x' for (t', r', 9' , 4>'), so that the requirement 
that the signal does not precede the turning on of the source requires G[x, x'\ =0 for t < t' . 

Note that it is not obvious that the appropriate spacetime dependence multiplying the (^-functions in 



the Green's function equation (31) is (rr') , and therefore we will justify this in appendix (TbJ below. 



When solving eq. (31) it is important to remember the following boundary condition. Since the 



Green's function is the field generated by a unit point mass, if we let ST describe a static point mass 
sitting at the origin, 



ST[x'} = SM5^ 3) [x] , SM/M < 1, 



(33) 



this merely amounts to shifting the mass of the central body by M — > M + 8M. Then we already know 



what to expect from the linear solution represented by the integral in eq. ( 32 ) . It should be the linear 



\\i-5M piece of the full IT[r] solution in eq. (14) upon the replacement M — > M + 5M. Remember that 



the mass dependence in the full solution of eq. (14) is contained entirely in r v via eq. (12). Perturbing 
M — > M + 5M in eq. (14) up to linear order in SM yields, 



n[r; M + SM] = n[r; M] + 5U[r], 



(34) 



8 



with 



8M i /rfijrf 11 



sum 



M p i 2irr v \ 6\/7r 



2F1 



117 r 3 
6' 2'6'~^! 



(35) 



Hence, when solving x'] below, we must obtain from eq. (32), 

dtG[x,x'] = SU[r], x = (t,x); x = (i',0). 



Mpi J-00 1 ' ' M p \ J_ 00 1 ^ ^ 

(The second equality follows from the time translation symmetry of the problem at hand.) For later 
use let us note that the small and large radius limits are, respectively, 

1 fr[§]r[|] 



dt G[x, x ] 



and (using eq. (95) below) 



27rr i; 



dt'Glx, x'] 



6\/vr 



r ->■ 



(37) 



4-7rr 



r — > 00. 



(38) 



III. RETARDED GALILEON GREEN'S FUNCTION IN BACKGROUND SOURCED BY A 
MASSIVE SPHERICALLY SYMMETRIC STATIC BODY 



In this section we solve, in different limits of physical interest, the Galileon retarded Green's function 



about the spherically symmetric background II [r] in equation (13). The first subsection summarizes 
all the results in a coherent manner, and in subsequent subsections we step through the derivation 
systematically. 



A. Overview of results 



Our solution of the retarded Galileon Green's function G[x, x'\ obeying equation (31) is composed 
of an integral over all angular frequencies uj and an infinite mode sum over all harmonics (£, m): 3 

/+00 j 00 +^ 

— e— Wu,J29l[ ur ,w / \ E (39) 
"°° 71 e=o m=-l 

We will make frequent use of the dimensionless variables 

(, = ur, = ur', £ v = ur v . (40) 



3 Strictly speaking, we need to specify a contour for the Fourier integral below, but since we do not need it in this paper, 
we shall leave this question for the future. 
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because we will shortly show that the radial Green's function depends on r, r' and r v solely through 



£, £' and £ v respectively. In eq. (39), the 17™ s are the usual spherical harmonics spanning a complete 
set of functions defined on a sphere of unit radius embedded in 3 spatial dimensions (the over bar means 
complex conjugation); they obey the eigenvalue equation 

L 2 Y™ = -£(£+ l)Y e m . (41) 

Because the background II is static, the Green's function reflects the time translation symmetry of 
the setup at hand. Moreover, spherical symmetry tells us the radial Green's functions gi do not depend 
on the azimuthal number m. 



The separation of variables method of mode expansion employed in eq. (39) reduces the problem 



of solving the linear partial differential equation for G[x, x'\ in eq. (31) to a linear second order 



ordinary differential equation (ODE) for the radial Green's function g#. Inserting the ansatz in eq. 



(39) into eq. (31), and using a Fourier representation of S[t — t'] and the completeness relation for the 
spherical harmonics, one may read off the ODE for g£ in frequency space by equating the coefficient of 
exp[— iu{t — t')]Y e m [8, 0]Yj, m [0 / , tfi'] on both sides of the Green's function equation. We then arrive at 



e 2 9| - \e,d, - ei + ^±^3) gt[t,?\ = (42) 



We have carried out a change of variables according to the rules in eq. (40); for the ei^,3 this amounts 
to simply replacing every r variable with its corresponding £ variable. That eq. (42) no longer depends 
explicitly on the radii nor on the angular frequency u) means that the solution for the radial Green's 
function gi, cannot depend on the lengths r,r',r v or frequency lv explicitly. 

The radial Green's function gi has a discontinuous first derivative at £ = £' because its second 
derivatives at £ = £' needs to yield <5[£ — £']/(££')• Therefore we need to distinguish between the two 
regions |£| > |£'| or |£| < |£'|. We therefore let r> and r< represent the larger and smaller of the two 
radii r and r', and define 

£>=wr>, £<=a;r<. (43) 
We now proceed to summarize the results for g£. 

Radiative Limit Of central importance in this paper, is the situation where the emitter lies deep 
inside the Vainshtein radius while the observer sits far outside (r 3> r v 3> r'). Taking the limit where 
one of the radii is much smaller than r v and the other much larger, more specifically r< <C r v and 
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rlr, 



FIG. 1: A Log-Log plot of the WKB "momenta" yj — U [r/r v ] (solid line) - see eq. ( [47] ) - as a function of the 
ratio r/r v . The asymptotics are: \/—U[0] = \/3/2 (long-dashed line) and y/—U[oo] = 1 (short-dashed line). The 
turning point, which is a global minimum, is at \/—U[r/r v = 1/2] = \JljZ. 



M>\ > l^| 3/2 , we obtain 



9t[W] 



(44) 



where J v is the Bessel function and the h^p is the spherical Hankel function of the first kind. 
In the high frequency regime \£ v \ S> £, the coefficients c'f^ are 



c, 



rad) 



3/4 



(45) 



,3/4 



where 



Xoo = jf (l - ^- di? ps 0.253 



(46) 



and 



U[#] = 



1 3(1 + 2t? 3 ) - 2 v W J (l + # 3 ) 



(47) 



4 1 + # 

and is plotted in Fig. @ below. The suppression factor of in ( |45[ ) indicates that high frequency 

Galileon signals are indeed Vainshtein screened, at least for small is. 
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In the low frequency regime, \£ v \ <C 1, the coefficients cl rad ^ are 



a 



(rad) 



(48) 



r[-f-|lr[Ml 



1 + 



cos[|(2£+l)7r] 
sinfrf] 



> 



2 2«+l 31 + gr [| + 3j r[ _ £] 

It is worth pointing out that although sin[7r£] appears in the denominator, and this expression 
contains T-functions whose arguments appear they could be negative integers; these I > 1 terms are 



all in fact non-singular. (This remark also applies to the related equation eq. (67) below.) We list the 
first ten C» rad ^s here, to 3 significant figures: 



t 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




0.324 


0.127 


0.0386 


0.0100 


0.00230 


0.000482 


0.0000932 


0.0000168 


2.86 xl0~ 6 


4.61 xlO" 7 



Notice from equation (48) that the £ = 0,1 modes do not contain r v ; this is a direct consequence of 
the fact that the leading order terms of the static Green's function in this same r> 3> r v S> r < limit 
(eq. (|67|)) do not contain any r v for £ = 0, 1. In fact, we may further compare the cubic Galileon radial 



Green's function g% with its non-interacting massless cousin 



-(Flat) 



ihf^^jdU] (see eq. ((77} and 



(78) below). Let us consider the non-relativistic limit <C 1, where the reciprocal of the characteristic 
frequency of the motion is much smaller than the characteristic distance of the source to the central 
mass M. In this limit, we may replace the Bessel functions with their small argument limits, and find 



~ -(Flat) 



ihP[Z>], 



91 = 9i 



(Flat) 



(49) 



and for £ > 2, 



A r > > r v > r<] 



9e 



(Flat) 



[r> > r v » r<] 



4d 



' V 



(50) 



fci W K>]- 



Here, K£ and k'^ are constants that depend solely on I. This teaches us that, while high frequency 
Galileon power loss is Vainshtein screened, low frequency signals generated from deep within the Vain- 
shtein radius of M are comparable to or even Vainshtein enhanced relative to the non-interacting 



massless scalar. The results of the radiative processes described in section (IV) will reflect these obser- 
vations. 
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The WKB High Frequency Limit When |£|, 3> max[l, £], we may apply the WKB ap- 

proximation. Let us first define 



(51) 



where it is worth noting that \/—U[0] = y/3/2 and \J — XJ\co\ = 1. The leading order WKB solution is 



m,?} = m/c v ]u[e/^]m 3 +e v )e(e 3 +e v )y 



(52) 



C++ 

I - exp [»& ($> -$<)] + -f— exp [i£„ ($> + $<)] 



where 



(53) 



- e -(f-i: 



> 



When max[l,£] <C <C the radial Green's functions become 



1 



-'^'J = ^7 U eX P &> " £<N + %" eX P [* ^> + ^< " 2 ^ X -)1 



(54) 



whereas when £ <C |£|, <C |fu| they become 

2 / 1 



9o 



exp 



+ 



exp 



^ (*>+*<) 



(55) 



In this limit, observe that the r> and i — t' dependent portion of the combination e * in the 



mode expansion eq. ( 39 ) is 



exp 



.V5 



V3 



(t - 



(56) 



and if we imagine a source located much closer to M than the observer is (r' <C r <C r„) the Green's 
function tells us the observer will receive purely outgoing radial waves. It is not entirely clear from 
the outset that this would be the case, particularly viewed from the curved spacetime picture (which 



we describe in section ( III C ) below) , because one may think that the radiation from the source could 
backscatter off the spacetime geometry and return to the observer, thereby mimicking an ingoing radial 
wave. (This scenario may in fact occur in the low frequency limit, where the longer wavelength of the 
Galileon waves may grow more sensitive to the curvature of the background effective geometry.) The 
phase r — (2/y/3){t — t') also indicates these outgoing waves, if they are propagating only in the radial 
direction, are superluminal because 2/y3 > 1. 
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Of particular importance is the WKB radiative limit, r> S> r v 3> r<. Here the radial Green's 
function is 



V2i 



V2i 



COS 



COS 



v/3 



+^ 



5-2* 







> 



(57) 



Large Mode Number It is useful to note that the asymptotic expansion of the Bessel function 
J u [z] exhibits an exponential suppression at large order v, for z < v. Therefore, the presence of 



Ji 



(i/4)(2£+i)[v / 3^</2] in eq. (44) means, at least for the radiation problem (r> 3> r v 3> r<) - the main 
object of this paper - we may neglect mode numbers much larger than |£<|, 

The Static Limit We may also obtain the zero frequency (static) limit of the Green's function, 
defined as 

-+oo 



G (static)[ fj ^ 



G [x = (t, x),x' = (t' , x')] dt 



(58) 



(It does not actually matter, because of time translation symmetry of the situation, whether we integrate 



with respect to t or t' in eq. (58).) As the name suggests, the static Green's function does not depend 
on time. We may interpret 

G (static) 

[x, x'] as the Galileon potential between two static point sources, 



both of unit mass, in the background U[r}. Putting the mode expansion in eq. (39) into the integral 



over all time in eq. (58) yields the mode expansion for the static Green's function 



oo +e 

G^[x,x>] = limu J Y J MZ,Z'} E ^[MW^'l 

W £=0 m=~l 

The exact result can be expressed in terms of hypergeometric functions as 



(59) 



^(static) 



+ 



2nr. 
1 




1 1 7 



EE 

t=\ m =- 



6y^ 

Yp[e,<f>]7p[9 ',4f] (r 
2£+l 

1 



6' 2' 6' 



2^1 



1 



1 t 7 I 

+ zz\ z: + 



6 6'2 2'6 3' 



2 Ft 



+ 



+ 



6 3' 2'6 3' 

eir [4(21 + 1)] 

Vvrr [1(21+1)] 



(60) 



r> \ 2 



2-^1 



1 



11 £7 £ r? 



6 6 2 2'6 



Via the identity in eq. (95) below, the £ = mode (the first line on the right hand side of eq. (60)) is 
equivalent to 



1 



47rr- 



114 

3' 2' 3' 



(61) 
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Also, since 



2 F X [a, P; T , z] = F[a tf^T^ (1 ~ Z V'^ [7 - a, 7 - ft 7 ~ « ~ I + 1; 1 ~ *] 



T[a]T[ 



r[ 7 _a-/3]r[ 7 ] 



2 Fi[a,/3;a + /3-7 + 1;1 - z] 



(62) 



r[ 7 -a]r[ 7 -/3] 

and using that l/r[— m] = if m is a positive integer or zero, one of the static mode functions can be 
written as, for t even, 



r J 



+ 



1 t 5 t r 3 



6 3' 2' 6 3' rl 



r 



r[ 



5-2fl 



r[-¥]r[ 



and, for £ odd, 



2-^1 



£ 1 
- + -. 

6 3 



„3l 



2' 6 3' 



r[ 



2-^1 



5-2^1 



f < + 2 1 r 3 
2'^ ; 2 ; + r! 



, (63) 



- 2 y 1+ r 3v^(;) 2 r[ _| ]r ^ 



2-^1 



1 ^ + 5 3 



6 



1 + 



(64) 



(Note that 2 F 1 [-£/2, . . . ; 1 + (r/r v ) 3 ] and 2 -Fi[-(^ - l)/2, . . . ; 1 + (r/r\,) 3 ] are, respectively, (^/2)th 
(even £) and (l/2)(£ — l)th (odd •£) order polynomials in the variable 1 + (r/r v ) 3 .) 

Taking the limit r, r' 3> r v hands us the Green's function to the Laplacian in Euclidean 3-space 



^f(static) ^ 



X, xf] 



Air\x — x'\ ' 



(65) 



plus corrections that begin at relative order (r„/r) 3 and (r v /r') 3 . This is to be expected, since far 
outside the Vainshtein radius, the central mass becomes irrelevant and we ought to recover the theory 
of a massless scalar in flat spacetime. 

Next, taking the limit r, r' <C r v leads us to 



^f(static) ^ 



X, xf] 



1 



yjrr' Jr% 



2nr v 



s/r/r v x — \J r'/r v x' 



+ 



mm 

6V* 



(66) 



a result obtained independently in [28] , and which will play a key role in our field theory based analysis 
of the conservative portion of the Galileon two body problem taking place in the background field 
II of M. Here x = x[0,(j>] and x' = x'[6',(j>'] are the unit radial vectors of the observer and source, 
respectively, and vertical bars denote the Euclidean length. We have expressed every occurrence of the 



two radii in eq. (66) as a small ratio, r/r v or r'/r v , to highlight the Vainshtein mechanism at work. 



In section (IIIC) below, where we shall view the Galileon propagating on the background II as a 



minimally coupled massless scalar propagating in a particular curved spacetime, we shall re-derive eq. 



(66) in an alternate manner. 
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For r> 3> r v 3> r< , we obtain 



1 1 

47rr> r> 




EE 

1=1 m =-e 



Yp[9,<i>]Yp[e',<!/] (y/rjr 



2£+l 



2r- 



(67) 



r[-f-^]r[|-|] / cos^ + iH \ 
v^r[-£] \ sinfyr^] ^ y ' 

plus corrections that begin at order (r</r„) 3 and (r t) /r>) 3 relative to these displayed terms. Notice for 
£ = 0, 1, r v drops out of these leading order terms. This tells us for the monopole and dipole terms, 
when the wavelength of Galileon signals are much longer than the Vainshtein radius r v , the Vainshtein 
mechanism becomes less effective. 

Well Inside Vainshtein Let both the observer and emitter lie well inside the Vainshtein radius, 
r,r' <C r v . In the low frequency limit, \£ v \ <C 1, the radial Green's function becomes, 



£„J_i V3£>/2 -V2Ji V3£>/2 



(68) 



for i = 0, and 



9d^']=i^J li2e+1) [V3U/2] (M^ +1) [VS^/2 



+ 



4^ +1 r [-I-|]r[| + |] 2 ^f 



i(-y - 1 <+§ 



3 f + | 7r 3/2 r 



2(1+2) 



J 



|(2M-1) 



V3£>/2 



(69) 



for £> 1. In the WKB limit, i.e. 3> max[l,^], the radial Green's function instead reads 

' i^e^M 1 ) [^>/2] J_i [V5e</2] £ = 



(70) 



M^e/Z Outside Vainshtein When both observer and emitter lie well outside the Vainshtein radius, 
r, r' 3> r v , we recover the theory of a minimally coupled massless scalar in flat spacetime, with the radial 
Green's function 



H %l + i) [v^>/2] J\(2i + i) [V3C</2] * > 



where in the low frequency limit \£ v \ <C 1, 



(71) 



a 



(hh) 



<0[£ v ] 



£ = 



^+i r [-f-|] r [-|-g](^[g( 2 ^+ 7r )]+ cot [g( 2 ^)] CSC M) p^ n 

2^+33 v / ^((2£-l)!!) 2 r[-^] 1 ^ U 



(72) 
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and in the high frequency limit 3> max[l,£], 



C, 



(hh) 



(73) 



^ _1 Ji_ e i'(HI)-MS»««) £> 



B. Solving The Radial Green's Function 

In this section we derive the results presented in the preceding section. In appendix (|B]), we review 
the relevant facts about solving Green's functions for linear second order ODEs, and also justify the 
(rr') measure on the right hand side of eq. (31). The algorithm for obtaining ge is as follows. 

General Solution of Radial Green's Function We need to first solve for the two linearly inde- 



pendent homogeneous solutions and Rf'iC] to the ODE in eq. (142]), namely 



(2)r 



e 3 <% - ei + 



£(£ + 1) 



e 3 )i?f' 2) [e] = o. 



(74) 



(1 2) 

In the notation of eq. (Bl), p2 = —e2 and p\ = —2es/C- Next we normalize the solutions R\ ' such 
that they satisfy 



i 



2 ' 



(75) 



Then the general solution to our radial Green's function is 



= c e R?>[z>]R? ) [e<] - (i - c t )R^[u]R?[^>] + c^R^m^ia + q 22 ^ 2 } K]^ 2) r], (76) 

where the constants Q, C^ 11 and C| 2 do not depend on £ nor but depend on The Q and 1 — Ci 
terms have discontinuous first derivatives and hence contribute to the coefficient of the <5-functions on 



the right hand side of the Green's function equation in eq. (42) 



For £ > 1, retarded boundary conditions and the demand for non-singular solutions will fix these 
constants uniquely. That is, we shall require that, whenever the observer is very far away from M, 
r 3> r v , and the source is closer to the central mass than the observer, r > r', then the observer ought to 
receive purely outgoing Galileon waves. Furthermore, on physical grounds, we will admit only solutions 
that are nonsingular when either the observer or the source is situated close to the central body. 4 For 



4 The observer placed close to the central mass will experience a Galileon force (oc 1/^/r) due to M that blows up as 
r — > 0; but here we are requiring that, as long as the observer is not sitting on top of the secondary source, i.e. the ST, 



she should not measure Galileon forces due to ST that grow without bound. 
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I = 0, in addition to regularity and the retarded condition, we shall also need to invoke Gauss' law 
applied to the curved spacetime Helmholtz equation (see eq. ( 131[ ) ) to fix these constants uniquely. 
In the following subsections, we will first solve for gg in the zero frequency (static) and high frequency 

(1 2) 

(WKB) limits. We will shortly also derive the R\ s in terms of Bessel and Hankel functions in the 
limits r, r' 3> r v and r, r' <C r v . This means we can fix the form of 7j£ in the limits r, r' <C r v , r, r' ^> r v 
and r> 3> r v 3> r< up the ^-dependent constants Cg, C\ x and Cf 2 . We will then proceed to fix these 
constants - at least within the low and high frequency limits, \^ v \ <C 1 and 3> max[l,£], respectively 
- by ensuring they agree with the static and WKB gg results in the same limits. 



Well Outside Vainshtein That we have just shown that the radial wave equation in eq. (74) 
reduces to that in Minkowski spacetime when r ^> r v implies that we may use the known solution there 
to read off the R\ ' s. The solution in flat (and, importantly, empty) Minkowski spacetime is textbook 
material, and we may represent it as 

5[t-t'-\x-x'\\ _ r du c _ luj(t _ t/) e^-*'\ _ 
Air\x — x'\ J 2tt Att\x — x'\ 

Then, using 



n = iw EE y "[ ' W, <f>'m<]hP[t>], (78) 



4tt\x — x , 

1 1 £=o m=-l 



we can deduce that the in the limit r>r„ must be a linear combination of jg and h^p 

The hp implement retf 
kind may be understood as 



The hp implement retarded boundary conditions, since the spherical Hankel function of the first 



fi) P , , I I d 



Note that only contains a factor of exp[+i£] and does not contain exp[— Thus, using the 

asymptotic expansion of the Hankel function for large argument, we see that 

e -Mt-t>) h W m H)W exp[-i,M-r)] (i + o[ (wrr i]) , (80) 

describes radially purely outgoing waves at unit speed propagating to infinity. A similar discussion 
shows that btp = (h^p)* implements advanced boundary conditions, and because jg can be expressed 

(1 2) 

as a linear combination of h\ ' , it describes a superposition of ingoing and outgoing waves. In the 



ansatz of eq. (76), we see that we have to choose r\ [£>] = hi [£>] and set Ci = 1 and C| 2 = to 



ensure retarded boundary conditions. (We are able to deduce from eq. (78) that h!jp and jg are already 



18 



appropriately normalized to obey the Wronskian condition in eq. (75) for r,r' S> r v .) This means we 



have determined g?[r, r' ^ r v ] to take the form in eq. (71 ). 



(1 2) 

Well Inside Vainshtein Let us now understand the forms of R f ' [£] evaluated close to the 



central source (r <C r v ). We exploit equations (16) and (17), keeping only the most dominant terms in 



ei2,3 (equations (28) through (30)), to reformulate eq. (74) as 



« + 2£ e 4 4£ 2 



)< ,2) [e] = o 



(81) 



One may rescale the solutions i?l 1,2 ' > [^"] = £ 1 / 4 ^ 1,2) [ v / 3^/2] and find that n { }' 2) [VS£/2] satisfies Bessel's 



equation 

where £ = \/3£/2- Noting that the Wronskian between J v and is 



(82) 



\N Hz) [J u ,HW] = J u [z](HWy[z] - (J u y{z]H(%] 



2i 



TTZ 



(83) 



and that 



Wr 



^J^^^iffW^l = v / iWr (z) \j v [z],HM[z] 



(84) 



we conclude that the two linearly independent solutions normalized to obey eq. ( 75 ) are 



4^ = 
4%] = 



2 ^3/2 ? 1(21+1) L 



2C 



^ 1 
3/2* J |(2f+1) 



%/2 



(85) 
(86) 



With these homogeneous solutions, the general solution for the radial Green's function, deep within 
the Vainshtein radius, r,r' Cr„, is 



1372 V« ^(M) 



2£ 



>/3e>/2 



j 



j 



i(2«+i) 



4 (2^+1) 



^</2]-(l-Q)iI« m) 



a/3^</2 



, ^(HH) rr(l) 
+ W ^1(2^+1) 



r(l) 



FT v 

i(2<+l) 



J i(2^+i) 



v^>/2 



(87) 
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Next, we recall the small argument limits (\z\ <C 1) of the Bessel and Hankel functions 

(z/2)i^ +1 ) 



J 



1(21+1)1 



r[\(2£ + 5)] 



:^+o[z 2 ]), 



TT 



+ 



(2£ + l) 



2 \ K 2 ^ 1 ) 



|(2^+l) (l+icot [f(2l+l)]) , 

r[i(2/ + 5)] ^ + ^JJ- 



The (2/z) u piece of the small argument behavior of H^\z\ = J v [z] + iiV^^] can be traced to J- U [z]. 

/TTTTN 

This implies that if we want a nonsingular solution as r < /r v — > 0, we must set Cg = 1 and C) = 
for £ > 1. For t = 0, however, both Ji/ 4 [\/3f/2] and £ X /4#W [^3^/2] are nonsingular in the small 

(2) 

radius limit; the former goes to zero and the latter to a constant. We must therefore write Rq [r< <C r v ] 
linear combination of these two functions. 
The Radiative Limit We may now fix the form of the radiative limit, r> 3> r v 3> r<, of This 
is crucial for studying the Galileon radiation seen by an asymptotic observer at r ^> r v generated by a 
source moving deep within the Vainshtein radius of the central mass (r' <Cr fl ). Our previous discussion 
leads us to the forms 



4 1} K>] ( ^ • v^i [V3^</2] + cf } • [V3e</2] 
[^[^Cf-^Ji^) [v^</2] 







£ > 



(90) 



1. Static Green's Function 



The static limit of the Green's function, as defined in eq. (58), requires a slightly different treatment 



from the time dependent case because it amounts to setting u = in frequency space, making the 



variables £, £' and in eq. (74) ill defined. We therefore work instead with the original radial variables 



r, r' and r v , in terms of which (42) becomes, in the uj — > limit 



e 2 cl + -e 3 (9 r 
r 



^+1) 



e3 )4 1,2|s) w = o. 



(Here e2,3 = e2,3[r, r v ] depend on r,r v , not £,£„.) We also express eq. (76) as 



=(s) 



r,r'] ee lim = Q< |s) [r>]i2? |s) [r<] - (1 - |s) [r<]i?f |s) [r>] 

til p(l|«) r„l d(1|«) r >] I /^22 p(2|s) 



(2k) r /l 



(91) 



(92) 
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and eq. ( 75 ) as 



e 2 [r,r v ] (< |s) [r](i?f |s) )'[r] - (fl£ W )'[r]J# |s) [r]) = 1 



(93) 



As we shall see, for ^ > 1, the radial static Green's function will be fixed once we demand that the 

~ (s) 

solutions are regular for all radii r, r . For 1 = 0, regularity is irrelevant; instead, g will be determined 



by ensuring that eq. (36) is obtained and, for r> — > oo, that the Green's function goes to zero. 
Multiplying both sides of eq. (91) by y/ r 3 (r 3 + r 3 ) (and dropping the labels) yields 



(r* + r$)K'[r] + 



4r 3 + r 3 (2 . £(£ + 1) 



Rr] 



-R[r 



(94) 



and this equation may be readily solved in Mathematica [30J. The general homogeneous solutions to 



the static radial mode equation eq. (91), normalized to satisfy the condition in eq. (93) are 



R^\r\ = 



{21 + l)r v \ r 



+ 



1 



6 ' 3' 2' 6 3' rl 



i+i 
r \ 2 



2F1 



1 I 1 £7 t_ r 6 
6~6'2 + 2'6 3'~rl 



_ (2£ + l)r v \r V/ 

Since 2-^1 [«, 7; z = 0] = 1, we see that the (1— Cf)R^ 8 ^ [r<]-R« [?*>.] term tends to \/r > /r v (r> / r < ) e / 2 
and the Cj 1 term in eq. (92) tends to (r 2 /(rr')Y^ 2 , as r/r v ,r' jr v — > 0. These two terms grow without 
bound, and therefore we must choose C( = 1 and Cl 1 = 0. 



Next, we use the identity 



2F1 [a, # 7; *] = ; [ - - 



T\0\T[<y-a 

mr[a-p] 
r[«]r[ 7 -/3] 



2-^1 



a, a + 1 — 7; a + 1 — /3; 



1 



(95) 



+ 



{-z)~P 2F1 



P, j3 + 1 - 7; /3 + 1 - a: 



to recast the product ■R^ 1 ' s ^[r>]12^ |a ''[r<] in eq. (92) in terms of 2-^1 [a, 7! — ( r t>/ r g) 3 ]- The two 
potentially divergent terms for £ > 1 are the ones proportional to 



? (2]s)r 





(■ 


r 2 , y 




r>\ 


l 

r v 


r <7 


r< 



7f] 



* 1 
3 + 6 



2C 22 T 



r 

^ 1 

3 + 6 



+ 2) 



+ 



2(1 + 2) 
3 



We may therefore collect the results 

Ci = 1, 



cf = 0, 



sin [|tt(2^ + 1) 

+ ^ ) { 

sin [^7r(2£ + 1)] J 

22 £\T[-±(2£+l)] 



(l + OK 2 /r|]) 
-1 (l + OK 2 /r|]). 



(96) 
(97) 



2^r[i(2^ + l)]' 



(98) 
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where we have used the r-function identity T[z]r[l 



7rcosec \7rz\ 



Notice that none of the regularity constraints apply for I = 0. In fact, Rq [r] is a constant. R t 
varies as \j2jr v \Jr jr v for small r/r v , and using eq. (95), we obtain the equivalent expression 



>(2|s)r 



(2[«)r 




mm 



TT 



2r 



2^1 



1 14 

3' 2' 3 : 



(99) 



which varies as \j2jr v (const. — r v /(2r)) for large r/r v . This implies that R [r] is regular for both 
large and small r/r v . To determine Co, Cq 1 and Cg 2 here, we recall the discussion towards the end 
of section (III]), that g^[r,r' = 0]/(4tt) must correspond to the coefficient of the 5M/M piece of 511 [r] 



in eq. (35); by spherical symmetry, the £ > 1 do not contribute to the solution generated by a point 
mass at the origin. This implies Cq = and Cq 1 = r [1/3] T [1/6]/ (6-y/F). When r> 3> r< » r v , by 
using the identity in eq. (95) on g$\r, r']/(47r), and setting the resulting to unity, we find that 
the only constant term (independent of both r and r') reads Cg 2 r 3 [l/6]/(9 • 2 2//3 v / 3~7rr„). Because we 



have already chosen the asymptotic boundary condition (see equations ( |14[ ) and (18)) that ip[r —> oo] 

ass 

see, at this point, that g^fr, r'] only depends on r> and not on r<. 



generated by a point mass located at some finite r' should approach zero, this implies Cq 2 = 0. We 



We may also arrive at the same result for gy\r, r']/(47r), without invoking the background solution 



IT, if we refer to the curved spacetime picture described in section (IIIC). In particular, the static limit 



of eq. (131), gotten by setting ui — > 0, translates to 

lim yj rr^d r g^ [r, r' = 0] 



r-s>0 



1-Cn 



1. 



(100) 



This immediately implies Cq = 0. Taking the r> S> r < ^> r v limit tells us that, as r> — > oo, we are left 
with 



ii 
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2r 2 [|]r 2 [i] 

7rr„ 




(101) 



Since we require (p[r — > oo] 
r[i]r[|]A/5F and Cf = 0. 



for any finite radial location of the point mass, we must have Cg 



ii 



This completes the derivation of eq. ( 60 ) . 

With the exact solution to the static Green's function in hand, we may now take the limits r, r' <C r v , 
r, r' 3> r v , and r> 3> r v S> r<. In the small radii limit, r, r' <C r v , we set the in eq. (60) to unity, 
and drop the subleading (rr'/r 2 )^" 1 " 1 -*/ 2 term relative to the dominant (r < /r > Y^ 2 \Jr < jr v term, to obtain 



^((static) |^ 



X, X 




e=o m =-e 



2£+l v^> V r > 



(102) 
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The infinite mode sum in eq. ( 102 ) may be collapsed into a closed expression by first summing over the 



azimuthal modes and then invoking the generating function of the Legendre polynomials. The result 



of the mode sum in eq. (102) is eq. (66). 



In the same vein, the large radii r, r' 3> r v static Green's function may be summed into a closed 
form by similar means. In fact, our solution is consistent with our earlier observation that Galileon 
dynamics reduce to that of a minimally coupled massless scalar field in Minkowski; for since we know 



its static limit is eq. (65), our Galileon static Green's function ought to reduce to the same in this 



asymptotic limit to lowest order in r v /r<. In detail, if one begins from eq. (60), applies the identity in 



eq. (95) to the and then proceeds to set the transformed 2-^is to unity - because their arguments 
will go as —(r v /r) 3 , which is very small at large radii - one finds a subleading term proportional to 
/(rr'Y +1 and a dominant term proportional to (l/r > )(r < /r > ) i . Keeping only the dominant term 



„2£+2 



and again converting the sum over spherical harmonics into one over Legendre polynomials, followed 



by applying the latter 's generating function, we reach eq. (65). 



As for the case r> 3> r v 3> r<, the result in eq. (|67|) can be obtained by starting with the exact 
ution in eq. 
2-Fis to unity. 



solution in eq. (60), but only applying the identity in eq. (95) to i?P' s ^[r>], followed by setting the 



2. The WKB Green's Function 



Next we consider the high frequency limit, |£|, \£ v \ 3> max[l,^]. We first rescale the mode functions 



via 



(£(e 3 +e)) 1/4 ' 



(103) 



so that the Wronskian condition in eq. d75b becomes [£](^f - C^)'^?^] = 1, and @ 



reads 



(104) 

Here U has already been defined in eq. (|47j), and we define 

t 1 _ 16£ 6 l(l + 1) + 4g 3 ^ 3 (5l(l + 1) + 6) + gg(4l(l + 1) - 3) 

16^(^3 + ^3)2 • i 1U5 ) 

For large |£|, ^ max[l,£], we see that the denominator of K scales as 8 powers of the large quantity 
1/5 ~ |£|, while its numerator contains six powers of 1/6 times terms of order unity, order I and 
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order £ 2 . This means that the largest possible scaling of K is that it goes as (£S) 2 <C 1. For £ = 0, K 
scales as 5 2 . Since U is of order unity, in the high frequency limit we may therefore discard K relative 
to U. 

Suppressing the irrelevant indices, we now seek to solve 

= -e 2 K"lZ] + UiaUnZ} ■ (106) 

(Here and below, we are introducing a fictitious parameter e that will be set to unity once the solutions 
to 1Z are obtained.) Observe that, viewed as a function of £, the U is very flat by assumption, because 
d$U[£/£ v ] = U'[£/£ v ]/£ v <C 1. This calls for the WKB method of solution, in which one uses the 

derivatives of U with respect to £ as an expansion parameter. We therefore pose the ansatz 

e (i/e)S[Q 

y/-U[£/€v] 



(107) 



It is important to note from eq. (47) that y/—U[£/£ v ] has no real zeros, though it has a singularity 
a t £,/£,v = — 1 and a global minimum at y/— U[£/£ v = 1/2] = y/2/3. Both in eq. (106) and the ansatz 
of eq. ( |107 ), e will turn out to count derivatives, so that 1/e implies an integral. Inserting eq. (107) 
into eq. (106) and setting the coefficient of each distinct power e~ to zero, the £ = term yields a 
relationship between U and 5', which we may integrate to obtain two solutions 



(108) 

The £ = 1 term gives a differential relationship between tvq)/ y/—U[£/£ v ], its first derivative with respect 

(109) 



to £, and 5' and 5". Through eq. (108), this gives 

T (o) 



constant. 



By setting to zero the coefficients of e e , for £ > 2, we find a recursion relation obeyed by rm, 



d^" 



(110) 



where — (or +) is chosen if we chose the + (or — ) sign in eq. ( 108| ). As advertised earlier, we see that 
every higher order in e contains an additional derivative with respect to £; and the 1/e in the phase of 



eq. (107) is the integral in eq. (108). 

For our purposes, we shall work only to lowest order in the WKB approximation, just involving S 



d 2) 

and T(o) ; the solutions R\ ' normalized to obey the Wronskian condition in eq. ( 75 ) are 



m exp 



3^1/4 • 



(111) 
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Let us pause to understand the large (r/r v 3> 1) and small (r/r v <C 1) radius limits. Examining Fig.Q 
reminds us that \/—U[r/r v ] is basically flat for large r/r v 3> 1 . Together with the limit \J —U[+oo] = 1 , 
we may infer, for some £o and £ such that > £o/£d ^> L 



^ / / 



i7[i?]dt? » & / s y^cTpjd^ + £ - e = c - & / (i - V-w) ( 112 ) 



We will justify below that we may now further approximate this integral by extending the upper limit 
of integration £o / (,v to infinity, 

rt/tv 



rvtv 

JO 



(113) 



where was defined in eq. (46). Similarly, for r/r„ <C 1, by the flatness of the potential y/—U[&] 
near i? = 0, we have 



Z/£v I .Pi 



We thus have 



cxp|±i^£ 



3->l/4' 



exp[±i(g-g„Zoo)] 



e/& « i- 



r/r v -C 1 



r/r„ » 1 



(114) 



(115) 



The r,r' <C r„, r, r' <C and r> 3> r„ 3> r< limits reported in equations (54), (55), and (57) follow 



from equations (52) and (115) once is computed. 

It is important to observe that, in this high frequency limit we are working in, i?^ [£] (the + sign 
solution) in eq. (Ill) is proportional to ^ 1//4 i^^ + i)/4 ["V^^ / 2] in the r <C r v regime and to hg [£] in 

the r > r„ regime. To validate this assertion, we merely need to compare the expressions in eq. ( |115 ) 

(2) (2) 
against the high frequency limit of the Hankel functions. (Likewise, R\ [£] is proportional to h\ [£] 

and £ 1/4 -^ ( ( 2 +1)/4 [v / 3£/2] in the limits r S> r v and r <C r\, respectively.) Retarded boundary conditions 



mean, therefore, that in eq. (|76l) we need to set Ci = 1 and Cf 2 = for all I > 0. At this point, our 



WKB radial Green's function solution takes the form in eq. (52). 

For £ > 1, may be fixed by regularity, demanding that the limit r < /r v — > yields a radial 

Green's function that is proportional to the high frequency behavior of t^J^ J(2^+i)/4[\ // 3£</2]. From the 



form in eq (52) and the large argument limit of the Bessel function this translates to the consistency 
condition 



7T 



cos 



tt2£ + 3 
2 4 



(116) 
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where \ is a constant. By converting the cosine into exponentials and equating the coefficients of 
exp[±i"v/3£</2] on both sides, this lets us solve for both and \. 

For t = 0, the radial Green's function can now be proportional to a linear combination of the 
high frequency limits of £< 4 J 1/4 K/3£</2] and £< /4 #{^[V3£</2], because as already discussed, both 



are non-singular in the zero radius limit. Let us consider setting r < /r v = 0, eliminating the £ 1//4 Ji 



/4 



term. If we now also take r> <C r v , and if we remind ourselves of eq. (115) and the large argument 



limit of H { J-\ we see that the WKB Cq term must match onto the high frequency limit of the 
yf&H^l [\/3f/2] term in eq. ([87]), i.e. C£ + oc C^ m when |&| > max[l,4 As we will 
discuss in section ( |III C ) below, the i = radial Green's function, which obeys an inhomogeneous 



Helmholtz equation in curved spacetime, must obey the Gauss' law in eq. (131): 



hm^c^o[^' = 0] 



1, 



(117) 



where we have taken the small radius limit of the effective metric in eq. (128). Because the 
Ji/4[\/3£/2] vanishes when evaluated at £ = 0, this condition applied to eq. (87) allows us to 
solve for 
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m p(HH) 
3/2 °0 



(1 + Q7T 

4d /2 ' 



(118) 



which in turn yields the Cq in eq. (53) 



In the preceding discussion, we invoked the flatness of W—U[d] at both large and small i9. Let us 
pause to quantify this flatness by examining the following expressions 

V3 



(119) 



and 



(120) 



The reason for the in front of the integrals is that, because these expressions occur in phases (e.g. 
exp[i£„3>g]), it is not sufficient for the integrals themselves to be much less than unity. We need the 
entire expression to be small, so that very little oscillation occurs. 

Computing the power series of Fold] and Too [l/<5] about 5 = allows us to understand precisely 



how good an approximation the expressions in eq. (115) are. For r/r v < 1 we obtain 



(121) 
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and, for r/r v 3> 1, 



Foo[r/r v ] = e,^^ + O [Utv/O 5 ] ■ 



(122) 



Therefore, replacing ro/r v with +00 in eq. (112) results in an error in the phase that is of order 



^(^/( wr o)) 2 "C 1, as long as we assume |£>| 3> ojtq 3> |£u| 3 ' 2 (i-e. the observer is far enough outside 
Vainshtein). Likewise, the approximation in for r < /r v <C 1, replacing \f—U with y/3/2, makes 

an error in the phase that is of order £ v (£</£v) 5 ^ 2 <C 1 as long as £ <C |^| 3 ^ 5 . 

3. Low and High Frequency Limits 

Because we were able to solve for the relevant mode functions (up to their overall ^-dependent 
normalization constants) in the regions very close to and very far away from the the central mass, 



we were able to determine the form of the radial Green's function in the limits r, r' 3> r v (eq. ( 71 ) ) , 



r,r' <C r v (eq. (87)), and r> 3> r v S> r< (eq. (p0[)). We may calculate the ^-dependent constants 
appearing in these expressions at least in the low frequency <C 1 and high frequency |£u| S> max[l,£] 
limits by matching them onto the static and WKB results we have obtained in the previous two sections. 

In the low frequency limit, we may simultaneously take the limit \£ v \ <C 1 and replace the Bessel 
and Hankel functions, and their spherical versions, with their small argument limits. The resulting 



expressions in equations (90), (87), and (71) can then be compared against the respective expressions 



in equations (67), (the small r/r v , r'/r v limit of) (|60[), and (65). 



Similarly, in the high frequency limit, we may simultaneously take the limit 3> 1 and the resulting 



expressions in (90), (87), and (71) can then be compared with the respective expressions in equations 



(57), (55) and (54). 



Notice that, upon these comparisons, for the radial Green's function evaluated deep inside Vain- 



shtein, eq. (87), Co = 1 in the high frequency limit, while it goes to zero in the low frequency limit. 



Finally, to arrive at equations (44), (68), and (70), it is useful to write H^X = J1/4 + ^Vi/4> and to 



use 



NJz] 



J u [z] COs[7T^] — J- U [z] 



Sin 7TZ/M 



(123) 



One may wonder why the result in eq. ( 71 ) is not simply the usual answer in flat Minkowski 



spacetime, i.e. Cp hfl ' = 0, since if both observer and emitter are very far from the central mass M the 
propagation of signals would not be expected to feel the presence of the central body M. The physical 
reason is that signals with wavelengths much longer than that of Vainshtein radius \£ v \ <C 1 indeed 
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cannot resolve r v very well - cf 1 ^ is proportional to some integer power of the small quantity - but 
once the wavelengths become much shorter than r v Galileon signals can resolve the Vainshtein scale 
and Cf hh ^ is not small but becomes a mere phase. 5 



We have now completed the derivation of the results in section ( III A ) . 



C. The Curved Spacetime Picture 

Before we move on to investigate the radiation seen by an asymptotic observer generated by the 
motion of matter close to the central source M, we would like to discuss an alternative perspective for 
the Galileon propagating on the background TL[r]. We shall also discuss the existence of a Gauss' law 
for the Helmholtz equation obeyed by the £ = mode of the radial Green's function. 



By a direct calculation, it is possible to view the Green's function equation (31) as that for a 



massless scalar wave equation in a curved spacetime. Specifically - recalling equations (28) through 



(30) - dividing both sides of (31) by y/e\ei&z, we arrive at 



a x G[x, x') = n x ,G[x, x'} = 6 ™}* a/ ? ] = VW/^WJ^W ~ f ] <^ (2) [x -#] ( 1 2 1 ) 

with 



100' 1 1/4 " rr ' 



_ 8# 3 / 1 + 1? 3 

m = TTi^y 3(1 + 2^) _ 2^mTTWy (125) 



and 



5 {2) [x-x') = <J[cos0- cos 0'](5[<£ -<£'], (126) 
= det g a/3 [x] , g' = det g Q/3 [x'\ . (127) 

We have denoted O x = Q a @[x]V x a V x p and n x i = Q a/3 [x']V x ia V x fp to be the minimally coupled massless 
scalar wave operator in a curved spacetime geometry given by the metric 

Q af3 dx a dxP = e^e|e 3 dt 2 - e\ e~ 2 K 3 dr 2 - ejelr 2 n AB dx A dx B , x A = (0,0) (128) 



Observe that the $j[r' /r v ] (eq. ( |125 )) occurring in eq. ( pi is (8/\/3)(r/r„) 3 in the small radii limit and 



unity at large radii. This is the Vainshtein effect at work: the magnitude of the point mass sourcing the 



5 A simpler toy example is to consider the theory of a minimally coupled massless scalar, with a spherical perfect 
absorber with radius R® centered at the origin of the spatial coordinate system. The radial retarded Green's function 
takes the form ih\ [£>](i<[£<] — X' l f 1 '[C<])- Perfect absorber here means the scalar field observed on the surface 
of the sphere is identically zero, thereby imposing % ~ je[^Ro]/h^ [ojRo]. At low frequencies, \uRo\ <C 1, X — ^ 
i(uiR ) 2e+1 1 \{2i - 1)\\{21 + 1)!!). At high frequencies, \uRo\ > 1, x ^ (V 2 )(l - exp[i(ir£ - 2uRo)]). 
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Green's function grows weaker the closer it gets to the central mass M, but goes to unity far away from 
it. On the other hand, it is somewhat puzzling, in this curved spacetime picture, that a source located 
nearer and nearer to the spatial origin tends to zero strength; for instance, we have already remarked, 
towards the end of section (jll]), that a static point source located at the origin must contribute to the 
background II[r] solution via a shifting of the mass, M — > M + SM. 

To understand this we first re-express eq. ( |124 ) in accordance with our decomposition in eq. (39); 
this means that we drop the integration symbol J dw/(27r), and replace 5[t — t'] and <5( 2 )[x — x'] with, 
respectively, e - lw (*~*') an d the spherical harmonic completeness relation. Because placing the point 
source at the origin means we have a spherically symmetric problem, this means only the monopole 
i = term in eq. (124) is relevant. Keeping the discussion general for now, let us merely assume the 
metric reads 



Q^ u dx^dx u = goodt z + g Tr dr z + 5ABdx A dx B , 
and is time independent. Our frequency space curved spacetime equation is then 



lj (-o;VM<? 00 5o[£,0] +d r ( vW^Sb&O 



Sr] 



(129) 



(130) 



which allows us to integrate over an infinitesimally small neighborhood about r = to obtain the 
normalization condition: 



lim wi 



lsls rr 0rSo[£,O] = l. 



(131) 



This indicates that, even though the measure $d{r/r v ]/r 2 (see equations (124) and (125)) tends to zero 
as r — > 0, a point mass sitting at the spatial origin must nonethless produce a non-zero (spherically 



symmetric) Galileon field, for otherwise eq. (131) cannot be satisfied. 



That we are able to obtain an equivalent curved spacetime wave equation to eq. (31 ), is in fact one 



way to justify the "measure" l/(rr') multiplying the 5-functions on its right hand side. From general 
arguments due to Hadamard |33] - see |34j for a review on Green's functions in curved spacetime - 



we know that solutions exist for the massless scalar Green's function equation (124). Since equations 



(31) and (124) are equivalent, this means the (rr') _1 is the correct measure. Moreover in this curved 
spacetime picture of the Galileon Green's function, we know that, when the observer at x and the 
source at x' can be connected by a unique geodesic, the retarded Green's function consists of the sum 
of two terms, namely 

e[t-t'] 



G[x, x'} 



4tt 



5[&x,x']y/ 'Ax,x> + @Wx,x']Vx,x> ) ■ 



(132) 
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Here, a x x i is half the square of the geodesic distance from x' to x. The first term after the equality 



describes propagation of Galileon signals on the null cone of the geometry in eq. (128); A x x > is related 



to the evolution of the cross sectional area of null rays emanating from x' to x. V[x,x'], known as 
the tail term, describes Galileons traveling inside the future null cone of x' . It is the solution to the 
homogeneous wave equation O x V = O x >V = 0, obeying non-trivial boundary conditions on the null 
cone of x' . We see that solving the retarded Galileon Green's function provides us information not only 



about the causal structure of Galileon signals but also about the effective spacetime in eq. (128). 



Via equations (16) through (20): we recover flat spacetime for the region well outside the Vainshtein 
radius (r 3> r v ), 



Q a/ 3dx a dx^ « rj^dx^dx" 



(133) 



and for the region well within Vainshtein (r Cr B ) we have instead 

Q /3dx Q dar : 



3 

' V \ 2 



' dt 2 - ^dr 2 - 3r 2 n AB dx A dx B 



2^/3 



(134) 



The 3/4 in front of dr 2 indicates that, deep within the Vainshtein radius of the central object, Galileon 
waves propagating in the radial direction are superluminal, since dr/dt = 2/V3 > 1. (We have already 
noted this, within the context of the WKB results, right after eq. (|55[).) In the same vein, it is worth 
mentioning that, if some method can be found to evaluate the infinite mode sum of the Green's function 



result for r, r' <C r v described by equations ( 39 ) and ( 68 ) through ( 70 ) , we should obtain the Hadamard 



form in eq. ( 132 ) and this would give us a deeper insight into the superluminal properties of Galileon 



signals near the matter source. 
Note that defining 



p = y V3r/2, p v = y V3r v /2 



(135) 



transforms eq. (134) into 







a pdx a dxP 



1 ( Pv 



(dt 2 - (2p) 2 5 ij dp i dp>) 



2\/3 V P 

where the Cartesian components of the spatial coordinates are 

p l = p(sin 6 cos 4>, sin 6 sin (f), cos 6) . 



(136) 



(137) 



This small radius curved spacetime metric provides an alternate means of deriving the inhomogeneous 



portion of the static Green's function in eq. (66). To see this, first re-scale the static Green's function 
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as 



G (static)[- - 



pi 



(138) 



Next, compute the static analog of (124), namely 



n x G^ tic) [x,x'] = \Q[x]e[x']\- 1 / i 6W[x-x') 



(139) 



using the small radii metric in (136). Thus, g( statlc ) is the Green's function to the Laplacian in Euclidean 
3-space 

d d 



p dp 1 dpi ' 



namely 



A p - ff ( static ) [p, jf\ = -A^g^^ [p, jf] = 5^ [p - if] . 



Remembering the rescaling performed in (138) and requiring that G( static ) [x, x 
fixes the general solution to take the form 

V3 ( pff 



G (static)^- 



4vr/)3 



p 1 



Xo 



Xi (p + P) - X2pp' , 



(140) 
(141) 

G( static )[x',x] then 
(142) 



where Xo,i,2 are spacetime constants; these Xo,i,2 terms are homogeneous solutions, i.e. O x and O x , 
applied on them give zero. The xo,i may be fixed by placing the source at the spatial origin, p' = 
or p = 0, and making sure that eq. (37) is recovered. (The xi may also be determined using the 



condition derived in eq. (131). First replace r with p; remember Q( statlc ) [x, 0] = lim^^o ^Qoi^ 0]/(4tt); 



and a short calculation yields vj^l^' 



pp 



jI/V3. Altogether, 



Xi 



1.) Finally, \2 can be fixed by 



demanding that, for r> —> oo, the Green's function goes to zero. (We do not compute Xo,i,2 in detail, 
since we have already obtained the exact result in the previous sections.) 



In equations (A9) and (A12) we derive the minimally coupled scalar field generated by n point 



masses in a generic curved spacetime in terms of the Hadamard form in eq. ( 132 ). The reason for doing 



so is that the Galileon field ip generated by the n body dynamics is, in fact, related to eq. (A12) - the 
curved geometry in question is eq. (128). 6 One cannot help but wonder if our mode expansion results 



in section (III A) can be utilized, at least in some limits, to extract the various portions (y / 'A XtX i, a XjX >, 



etc.) of the Hadamard form in eq. (132). We leave these questions for possible future work. 



It is not exactly the same expression because, in curved spacetime, the element of proper time is ds = */ g^ v x fJ -x"dt, 
whereas in this paper we are working in flat spacetime and the element of proper time is ds = rj^ lv x> i x v At. 
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IV. RADIATIVE PROCESSES 



In the following two sections we wish to examine the Galileon waves produced by matter in motion, 
within two concrete scenarios. The first is the radiation made by acoustic waves propagating on the 
surface of the large central body of mass M with radius Rq. The second is the radiation created by 
the movement of compact bodies orbiting this central body; for instance, this could describe our solar 
system, with planets orbiting around the Sun or a highly asymmetric mass ratio binary star system 
capable of also generating gravitational waves that could be observed by upcoming gravitational wave 
detectors. 

We initiate the discussion by defining radiation to be the piece of the Galileon field that carries a 
non-trivial energy-momentum flux to infinity. Quantitatively, we only wish to consider the portion of 
ip that contributes a non-zero power (per unit solid angle) at infinity, 

d W 

— — ee lim r 2 T 0r = - lim r 2 d r ipd t ip. (143) 
(We have used eq. ([9]).) As we will see more explicitly below, both dfip and d r ip contain a power series 



in 1/r, beginning at 1/r. Even though the full Galileon field is II ee II + <p, according to eq. (19), the 
derivative of the static background Galileon field goes as II' [r — > oo] ~ 1 /r 2 , it does not contribute to 
the asymptotic power we are currently after and hence would be neglected in the following discussion. 

To highlight the importance of the non-linear self interactions of the Galileon field, we will also 
compute the radiation signature in the same setups arising from the minimally coupled massless scalar, 
so that we can compare the prediction of the two theories. That is, we will also consider the theory 
described by the following action 

^/ A (l ( ^ +s 3^/* (i44) 

Compare this with the Galileon action Sn + Sm + 5S in equations Q, ([7]), and Q. There, we had 
to first solve II sourced by M and then proceed to perturb around it. Here, because the H-theory is 
linear, the complete field H is gotten by superposing the field spawned by each individual source. 

When we compute the radiation generated by the surface waves on M, we will assume that these 
waves are driven by external (non-gravitational) forces, so that M itself can be considered to be a static 
source of the scalar field(s) (falling off as 1/r 2 ), and do not produce any radiation due to backreaction. 
The only source of radiation there is 5T describing the waves themselves. When we compute the 
radiation of the n-body system, however, we need to remember that, even though M ^> m a , and hence 
the location of M never deviates far from the center of energy of the system - the entire system is 
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held together by (largely) conservative gravitational forces. In this case, we shall see it is important to 
include M in the radiation calculation so as to enforce the conservation of linear momentum. 



The radiative limits of the Green's function for ip (equations (44) through (148])) and for H (equations 



(|77j) and (|78j)) both take the form 

dw 



G[x,x'} 



-iu(t-t') 



2vr 



>J > 



(145) 



where we have chosen to write the orthonormal angular mode functions as a linear combination of the 
usual spherical harmonics 



r>m _ Tin v" 
U t - ( £ )L n ,Y £ 



(146) 



(The * in (f2™)* represents complex conjugation.) Here ^L m n , is a unitary (2^ + 1) x (2^+1) matrix, so 
{0™} is as good an orthonormal basis as the spherical harmonics. For the minimally coupled massless 



scalar S, equations (77) and (78) tell us 



while eq. (44) gives, for the Galileon cp, 



n 



«\<p) 



uC. 



(rad) 



/ 3?</2 



(147) 



(148) 



where ao = —1 and og = 1 when i > 1. If we decompose the matter source in the same way that we 
decomposed the Green's function, 



5T\x] 



-iu)t 



2vr 



(149) 



£,771 



then by the orthonormality of both the exponentials and the angular mode functions, we may translate 



eq. (32) into the following solution for the cp or H field evaluated at (t,r,9,<p), sourced by 5T/M P \ 



ip[x\,r,[x\ 



roc 

7Z Jo 



u, r 



(150) 



The asymptotic behavior of the spherical Hankel function then implies that the scalar field, for a fixed 
angular frequency uj, is indeed proportional to a finite power series in 1/r, and the series begins at 
1/r. This leading order 1/r piece of the time and radial derivatives is precisely what we call radiation, 



because when inserted into eq. (143), the factors of r cancel and what remains is a finite amount of 



energy transported to infinity. Any part of the fields containing a higher power than 1/r would make 
a contribution to r 2 T 0r that decays to zero at infinity. 
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Moreover, since taking a time derivative brings down a — iuj and taking a r-derivative - to lowest 
order in 1/r - brings down a we deduce that the radiative part of the fields obey the relationship, 



iptf) radiation — [dr'v) radiation 



with H respecting the same equation. Hence it suffices to display only the time derivatives, 



(151) 



d t <p, d t 



1 



pi 



duj 
2^ 



l + O 



dr'r u Kf[t;']pT[u,r']. (152) 



Because it contains an r-dependent exponential e^, one may worry that the Fourier integral in eq. ( |152 ) 
would somehow introduce additional terms that go as 1/r (after the w-integral has been performed), 



and therefore that dropping the higher powers of 1/r in eq. (152) at this stage is premature. However, 



we will see from the integral in eq. ( 154 ) below that the e*^ drops out of the expression for total energy, 



and therefore r does not take part in the resulting w-integral; in particular, the higher powers of 1/r 



that have been discarded in eq. (152) would indeed decay away once we take the r — > oo limit. 



In general, we expect the Fourier integral in eq. (152) to be extremely difficult to evaluate. We will 



instead consider what the spectrum of radiation emitted from a particular system is. The spectrum is 
indeed a physical observable, since observations usually take place over enough cycles of the radiation 
field for a Fourier analysis to be done. 

The total energy per solid angle radiated to infinity is the integral 



dE 

da 



lim 

r— >oo 



r 2 d r (fdt(fdt. 



(153) 



(The same expression holds for H.) A few standard Fourier identities allow us to express the total 



scalar energy radiated per unit solid angle, from eq. (152), as the following integral over all angular 
frequencies 

2 



dE[tp or S] 
dfl 



1 



Ml 



dco 
2^ 



e,m 



dr'r' 2 n<lt']pTlco,r'] 



(154) 



In the subsequent two sections, we perform the decomposition in eq. (149) for surface waves on the 



spherical mass M as well as that of n compact bodies orbiting it, and proceed to apply them to eq. 
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A. Surface Waves On Spherical Body 



In this section we will describe surface (acoustic) waves propagating on the large central mass M of 
radius Ro by 

M 

STq[x] = —5RS[t-Rq]. (155) 

v e 

where Vq is the volume of M, 

V & = ^Rl (156) 

Denote by R[t, x] the radius of the mass M at a given time t and direction (0, (ft) from the spatial center 
of the coordinate system. For a non-relativistic system, which we shall assume is the case in this section, 



the trace of the stress-energy tensor ST primarily describes its mass density STqq. Then eq. (155) may 
be interpreted as describing surface waves of very small fluctuations SR[t, x] = R[t, x] — Rq around the 
mean radius Ro, on an otherwise perfectly spherical body. We will decompose these undulations SR as 

oo I 

SR[t, x]=RoJ2Yj (Af[x\af cos [Sift + &?] + B?[x]bT cos [n' e m t + VP] ) , (157) 

£=l m=0 

with \af\, \b™\ <C 1. The "rotated" spherical harmonics {A™,BJ 1 '} are defined as, whenever 



Af = i m (Y e m + Y e - m )/V2, 

%T = i m+1 (Yr - Y e ~ m ) /y/2, (158) 

whereas for m = 0, 

1?^° = (!?)*, B^° = 0. (159) 

Because Y^ m = (— ) m Y^ m , the A™ and B™ are real. This is the primary reason for using them instead 
of the usual Y^ m , because now aj 1 and bj 1 can directly be read off as dimensionless amplitudes of a 
particular mode of vibration, with respective oscillation frequencies OT 1 , f^ m , and phases $™ and fy™. 
For technical convenience, we will assume all oscillation frequencies Q'™} are distinct and positive. 

It is also worthwhile to observe we have not allowed a monopole £ = term in the infinite sum eq. 
(157): in particular, f d 3 xST[t,x] = 0, and the total mass is a constant. While our scalar field theories 



do not enjoy any sort of symmetry giving rise to a conservation law for the associated scalar charge, 
the coupling to the scalar fields here is of (sub-)gravitational strength, and thus very weak. Therefore 
the requirement from the known laws of physics that mass is a conserved quantity, in the Minkowski 
spacetime we are working in, therefore takes precedence. 
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By a direct calculation, one can check that our surface waves have the following decomposition 



Yl E f^n^K (e _i *™5[u; - + e^S[oj + fiJT 1 ]) + B^[x]bf {e~^ 5[lo - fij 



£=1 m=0 



(160) 



Here, the analog of the p™s appearing within the formula in eq. (154), are the coefficients of 

. | 2 in eq. (154), we en- 



f(diu/(2Tr))e- iu}t A™ and J(du/(2ir))e~ iu;t Bf. When taking the square 



counter cross terms involving the <5-functions. However, since we have assumed that all frequencies 
are distinct and positive, the arguments of the <5-functions cannot be simultaneously zero unless the 
frequencies are in fact the same. This collapses the summations into a single I- and a single m-sum. It 
remains to deal with squares of the form ((2tt)5[lo ± w']) 2 , with u/ being one of the 0™s or We 
will treat one of the (2it)5[uj ± uj'\ as the total duration of time, since 



(2tt)5[u) = 0] 



dt lim e~ iujt 



total time elapsed. 



(161) 



Dividing both sides of eq. (154) by total time, i.e. (2tt)5[uj = 0], then gives us back power radiated in 
scalar waves per unit solid angle. 

Minimally Coupled Massless Scalar The result for S, from eq. (147), is then 



dm 

dtdfi 



9G N M 2 g £ 



nR 2 



\A?[x]a?n?Roj t [n?Ro 



+ 



(162) 



(We have exploited the fact that, because ji[z\ is z l times a power series in z 2 , \je[—z)\ = \je[z]\.) 

For non-relativistic systems, a substantial subset of the oscillation periods ~ l/0» m ought 

to be much longer than the light crossing time ~ Rq. Therefore we can treat the small dimensionless 
quantities £l 1 J l Ro,Q'J n Ro <C 1 as expansion parameters. This leads us, in this non-relativistic limit, to 



dm 

dtdQ 



9G N M 2 



EE 



ttR 2 ^ ((2£ + l)!!) 2 

1=1 m=0 u 'II 



Af[x]af (QfRo 



+ 



(163) 



Galileons For the Galileon ip, because the Vainshtein radius depends on A in eq. © - a fr cc 
parameter in this paper - it is not necessary that iVpr v , ^l'pr v C 1 < f. This prompts us to write 
the energy output dE^/dtdO, due to Galileons as a function of the mode numbers (£,m) instead. We 
will assume that the surface waves are non-relativistic, and thus the ratio of the radius to that time of 
oscillation is a small number (i.e. Q^Ro, Q'^Ro <C 1), so that we may replace the Bessel and Hankel 
functions with their small argument limits. 
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When Ufr v , Vl!™r v < 1 and £ > 1, 

d£ £ m M _ 288G N Af 2 Ro ( r [-f - i] T [| 



dtdn 



2 3 ^v [f + 1] r [| + 1] rH 



'cos [|vr(2£ + 1) 
sinf7r£l 



+ 1 



(164) 



Notice that for the t = 1 case, the Vainshtein radius r„ drops out; low frequency oscillations 
are thus unscreened to the lowest order. At higher than dipole order, £ > 2, by comparing the 



(Ro/r v )(nfRo) e+2 (nfr v Y in the Galileon result in eq. pil) to the analogous (flfRo)' M+2 in eq. 



\2£+2 



( fl63l ) for the non-interacting massless scalar, we find that Galileon power is actually enhanced (up to 
^-dependent numerical factors) by the ratio (r^/iZo)^ 1 1. 
When Qfr v , n' e m r v > max[l,4 and £ > 1, 



dE?[<p] _ G N M 2 
dtdn 



e , 9 

32+4 



4^r[f + 



f] 2 



A?[x)a? 



2 (Of Ro 



d+3 



+ 



Bnm 



(n' e m R Q 



(165) 



For a fixed I mode, we may take the ratio of the S results in eq. ( fl63l ) to the ones here for high 
frequency oscillations: up to numerical ^-dependent factors, we obtain (f2™i?o)^ _1 (0™r t ,) 3 / 2 . For 1 = 1, 
this ratio is (f2™7\,) 3 / 2 3> 1, telling us the dipole term in eq. (165) is Vainshtein screened. However, 
in this non-relativistic limit, Q^Rq <C 1, once £ is large enough that (O^o)^ 1 ^^) 372 < 1, we see 
that high multipole Galileon radiation becomes Vainshtein amplified relative to their non-interacting 
cousins. We also note that, while low frequency oscillations in eq. ( |164[ ) contain integer powers of 
frequency, here the radiation spectrum contain (Q^Ro) 3 ^ 2 and (Q'^Rq) 3 ^ 2 . 



The expressions in equations (164) and (165) came directly from the Bessel and Hankel mode func- 



tions satisfying the original wave equations of <p and S. Together with the presence of the combinations 
{l™Ro and Q'jT'Ro, these facts teach us that radiation generated by the surface vibrations of the massive 
object M directly probes not only the vibrations themselves but also the dynamics of our field theories 
- in the Galileon case, it carries information about the theory operating deep within the Vainshtein 
radius of M, where the self interactions of the full II theory are dominant. 

At small ratios of the Vainshtein radius to oscillation time scale (eq. ( |164 )), we see that the Galileon 



radiation spectrum, like its H cousin in eq. (163), contains only integer powers of the oscillation 



frequencies. However, for very large Vainshtein radius to oscillation time scale ratios (eq. ( 165| )), the 
power emitted begins to contain fractional powers of angular frequencies; this indicates there must 
be a change in the spectral index if one is able to probe Galileon radiation over a broad bandwidth. 
This can be traced to the Galileon radial mode functions </(i/4)(2^+i)[V / 3£/2] within the Vainshtein 



37 



radius, r <C r„. We will witness this phenomenon again in the following section on the n body radiative 
problem. 

B. n Point Masses Orbiting Within The Vainshtein Radius Of Large Central Mass 

In this section, we consider an arbitrary number of compact bodies of masses {m a \a = 1,2, ... ,n} 
orbiting around the central object M? We will assume this system is held together by only gravita- 
tional and scalar forces, and we will further approximate these compact bodies as point masses, with 
spatial position vectors {y a \a = 1,2,..., n}. If we work within the non-relativistic (i.e., slow motion) 
approximation, valid for a wide range of astrophysical dynamics, including that of our solar system - 
the Galileon and gravitational interactions are described by 

|>/-H(^)--)(^M- 



The effects of gravitation are encoded in the proper time element ds a = dt\f g^y^y", which we have 
expanded in powers of velocities and the graviton field h^ u /M p \. We have assumed gravity is weak, 
9[iu = Vfiv + h^v/Mpi, where {h^/Mpil <C 1. The virial theorem tells us that the potential h^ u /M p \ 
scales as the square of the typical velocities v\ = (dy a /dt) 2 , so that the "+..." in eq. (166) can be 



understood to scale as and higher. (The scalar potential S/M p i would scale similarly; but the static 
portion of ip/M p \ would be considerably weaker because of Vainshtein screening.) 

Because we are interested in Galileon and not in gravitational radiation, we may "integrate out" 
the gravitational field hoo/M p \. To the lowest order in the non-relativistic expansion, this amounts 
to replacing the hoo/M p \ evaluated on the ath point mass world line with the gravitational potential 
exerted on m a by M; the potentials due to the other compact bodies, as long as they are distant 
enough, scale as rrib/M, b ^ a, relative to that due to M, and hence are subdominant. (We ignore the 
possible self- force contribution.) The Newtonian energy per unit mass E a /m a = (1/2) (y^ + h^/M-S) is 
a constant (up to corrections of 0[w„]) - this means we may replace (l/2)(y^ — hoo/M p \) in eq. ( |166 ) 



with y 2 — E a /m a - At this point, what is relevant for the Galileon radiation problem is the interaction 



7 The gravitational waves analog to this section can be found in |31j . The conservative aspect of the n body (weak field) 
gravitational problem has a long history; see, for instance, |32j and the references within. 

8 The need for including the gravitational potentials, in order for the ensuing analysis to be consistent with energy 
conservation, has been emphasized in [29], and our discussion here overlaps with that treatment. 
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with the non-relativistic expansion denoted by 



v 3 v 2 



m 
M 



(168) 



By a direct calculation, the trace of the associated stress-energy in eq. (167) is 



5T m [t,x} = Y / ^a J dt'F a [t']S[t-x°[t']]S^[x-y a [t'} 



(169) 



We rewrite this as 



iM = J 



dtu 
— < 
2vr 



oo I 



-Ujt 



E E YrmpTi«,r], 



(170) 



t=0 m=-£ 

where, by a change of variables from proper time to coordinate time 



n" 1 

Pi 



[oj,r)=Y / m a J dt'Falt'jj' 



(171) 



r\ya[t']\ 

With this pj 1 , the total scalar energy emitted by n bodies orbiting around the parent body M is 



therefore given by (154). 
dE[<p or H 



327TGN 



dcj 
2tt 



a=l 



E^ / dt'i^v^'E E Yr[x]Yr[ut']](-iYn<[u;\y a [t'} 



=o m=-e 



Minimally Coupled Massless Scalar Using equations (147), (168) and 

oo t 



-ik-x 



47r E E (-i) e je[\k\r]Yp[k]Yrm 



(172) 



(173) 



i=o m =-e 



eq. ( 172 ) becomes 

dE[E] _ 2G N f du: 



a=l 



E m a / dt'e^'dt, {F a [t'] exp [-iuix ■ y a [t']] } 



(174) 



where we have also used the fact that every integer power of uj occurring within our integrand may be 
replaced with a time derivative (namely, idf) acting on t' dependent factors. Let us now Taylor expand 
exp [— iujx ■ y a [t']], and convert each additional power of u into an additional idy. Roughly speaking, 
each time derivative should scale as 



d t ~ w ~ v a /r a , 



(175) 



where r a is the typical orbital radii of the compact bodies in motion. 
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We have previously highlighted that the n-body system under consideration is held together by 
(largely) conservative forces, and it is therefore important to include the central body M in our radiative 
calculations to respect the conservation of linear momentum. This is an appropriate place to consider 
how the backreaction on the central mass M by the compact bodies orbiting it affects our analysis here 
and below. What we shall do is to extend the sum Yli< a <n ~~ ^ Eo<a<n to include that of M itself, via 
the definitions 

yo = spatial location of M, mo = M. (176) 

For the minimally coupled scalar case at hand, this is nothing but the principle of superposition, since 
the theory is linear. For the Galileon case, this amounts to including in the matter perturbations 



6T m in eq. (167) the time dependent multipole moments of the central mass M, induced by the 
gravitational and scalar forces of its n planetary companions pushing it away from the center of the 
spatial coordinate system. (That all our results in this section depend on at least two time derivatives 
of the spatial coordinates of M and the n compact bodies corroborates this interepretation.) This is 
subtly different from the non-self-interacting H case because, to capture the Vainshtein mechanism, we 
first had to assume that M was motionless so that we could solve for the background IT it generated, 
before proceeding to compute the Galileon Green's function. It is for this reason we have phrased our 
discussion of including M in the sum over a in terms of a backreaction on the motion of M. More 
quantitatively, the stress energy of the central mass is given by M5^[x — yo], and we may Taylor 
expand it about yo = 0. The lowest order term is M5^ [x] , which is the source of the background field 
IT; we are thus treating every higher term in the expansion (which is necessarily proportional to powers 
of y a [t]) as part of the matter perturbation 5T m . 

By Taylor expanding exp [— icux ■ y a [t']], followed by converting all cj's into time derivatives, we find 
that the non-relativistic scalar energy loss per unit solid angle from our n-body system is given by the 
expansion: 



dE[E] _ 2G N f doj 



dQ n / 2vr 



n 



2^m a J Me ix + 2 +U 

a=0 J \ 

The importance of including M in the sum over a is now manifest, for the first term involving the sum of 
all forces must yield zero, ^ a m a (d 2 y a /dt' 2 ) ■ x = 0; in a spacetime translation symmetric background, 
linear momentum is conserved and Newton's third law must be obeyed. Therefore S radiation really 
begins at O [vjM . 



Va,Va M 



(177) 



Observe the hierachy: |yb| <SC \y a \ r v , where a = 1,2,3, ... ,n. 
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As a check of the formalism here, in eq. (A20) below we shall re-derive the analog of eq. (177) 



directly from the position space Green's function in eq. (77), but (for simplicity) without including the 
gravitational potential in the proper time element. 10 We see that, at this order in the non-relativistic 



expansion, including the gravitational potential merely changes the — (l/2)dy^/dt' in eq. (A20) to 
-dyl/dt' in eq. (fir?). 



Before we proceed to the Galileon case, let us note that we could have obtained eq. (177) directly 



from the infinite £-sum in eq. (172), if we Taylor expand the spherical Bessel functions and use the 



explicit polynomial expressions for the P^s. (The reason for collapsing the infinite f-sum into an 



exponential in eq. (174) is to emphasize that, because exp [— iujx ■ y a [t']] admits a Taylor expansion 



in integer powers of co, the radiation spectrum in (174) depends on frequency solely through time 



derivatives acting on the y a s; this will not be the case for Galileons.) Because we are seeking an answer 
accurate up to Ofu 3 ], by counting powers of u), we may infer that up to the 1 = 2 terms of the sum 



are needed. Specifically, the m a d 2 y a /dt' 2 term comes from the leading order piece of the £ = 1 term; 
whereas the d 3 (x ■ y a ) 2 /di' 3 from the leading order piece of the I = 2 term; and the dy^/dt' comes 
from the £ = term with F a included. Furthermore, there is cancellation between the £ = first order 
correction term involving (o;|y a |) 2 and that from the £ = 2 term. 

Galileons In parallel with the treatment for 5, we will assume that our astrophysical n-body 
system is non-relativistic, so that the orbital time scale 1/oj is very small compared to the light crossing 
time ~ \y a \i allowing us to replace the Bessel and Hankel functions with their small argument limits. 

For low frequencies, \£ v \ <C 1, we have 

dE[<p] 2G N 



d^d(w/2vr) vr 



o=l 



£m« / dt'e^'Y^^iFalAMAt']} 



1=0 



(178) 



where 



M e [t'] 



ny a [t'}\) 2 



+ 



(179) 



\y*[t'l\ 



{r v \y a [t'W P,[x-y a [t'W 



A(2m) r[-f-!]r[§-§] 
2 3 ^r[f+i]r[f+f]rH] 



cos[|tt(2^+1)] 



+ 1 



> 



and we have converted the spherical harmonics sum to one over Legendre polynomials. The presence 
of dfi implies the magnitude of the £th channel is suppressed by \/r/r v (r/rY^ 2 relative to the lowest 
£ = mode, where r is the typical orbital radius and r is the typical orbital time scale; the factor 



10 Be aware that the frequency space analysis here does not capture the dependence on the approximate retarded time 
t' = t-r. 
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of (r^/r)^/ 2 further suppresses the power loss because the motion is highly non-relativistic. (We have 
already noted that the leading monopole and dipole terms in the radial Galileon Green's function in 
the non-relativistic limit matches that of its non-interacting cousin.) 

Developing the non-relativistic expansion up to requires the monopole, dipole and quadrupole 

terms (£ = 0, 1, 2): in the low frequency limit, we gather 

d 2 y a 



dE[<p] 



2Gn 



dftd(u;/(27r)) vr 
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dt'e iujt ' x 



dt' 2 



id 3 yg dgg vzn-pm d 3 

4 dt' 3 dt' 40^ 3 / 2 dt' 3 
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(180) 



Comparing equations (177) and (180) informs us that very low frequency Galileon waves travel essen- 



tially unscreened, even though they are generated deep within the Vainshtein radius of the system. The 
y/ r v/\ya\ factor in eq. ( 180 ) may even yield Vainshtein enhancement relative to its cousin S. Moreover, 



we shall now argue that unlike eq. (177), the first term ^2 a m a (d y a /dt ) • x here in eq. (180) and 



in eq. (183) below, is small - it scales as m/M - but is no longer exactly zero. The primary reason 
is that, while gravitational forces between any two objects obey Newton's third law, Galileon forces 
between the compact bodies (in the non-relativistic limit) do not. This is because the background II 
does not respect spatial translation symmetry; this statement can even be checked explicitly by taking 



the gradients of, say, the static Green's function in eq. (66) with respect to both x and x , and noting 



they do not give equal and opposite spatial vectors. It does turn out, however, that the force between 
each compact body and the central mass is equal and opposite: without loss of generality we may 
consider some small mass m lying on the positive z-axis. The problem is now cylindrically symmetric, 
which tells us the force on M due to to; and the force on to due to M, must both point along the 
z-axis. The J dt{M or m)Yl/M p \ coupling tells us the z-component of the force on to due to M, to 
leading order, is simply (m/M p i)Il' = —Mm/(2TrM^r^' 2 ^/r') (refer to II' in eq. (|16|)), where r' is the 
radial location of m. The force on M due to m is, in the non-relativistic limit, given by first invoking 



the cylindrical symmetry to replace \\/r/r v x — \Jr' jr v x'\ in eq. (66) with \ y/r/r v — \fr'/r v \, and then 
computing lim r ^o(Mm/M 2 l )d r G^ statlc ^[x, x'] = Mmj (2TiM^\j 2 \fr')} x Note that even though New- 



11 The reader concerned about the domain of validity of the ip solution generated by m, whose gradient is responsible for 
the force acting on M, can perform the following order-of-magnitude check. Replace the II in the Lagrangian density of 
the action in eq. Q with the total field of M and m, i.e. II -> 11+ (m/M p i)G (static) [x, x 1 ], with G (static) [x , f '] given by 
eq. (661. After expanding about r = 0, divide the dominant piece of the quadratic-in-G ( - statlc ' ) portion of the resulting 
Lagrangian density by the dominant piece of the cubic-in-G' statlc - ) term. For fixed r', the radial location of m, one should 
find the ratio to go as Mr'/(mr) 3> 1 - i.e. the linear solution offered by c( statlc ) should be an excellent description of 
the force of m on M. In the same spirit, one may also expand this same quadratic-to-cubic ratio about r — r' and find 
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ton's third law is not obeyed between the compact bodies, because the theory we started with in eq. 
([5]) was defined in flat spacetime, total linear and angular momentum must still be conserved. What 
must happen is that the radiation generated at 0[v%] carries away some of the linear momentum. 
At high frequencies but non-relativistic orbital speeds, i.e. r v S> l/|w| 3> |y a |, 

2 



dE[(p] 



dftd(u;/27r) 
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(181) 



+ ... 



(182) 



■ 2f+r, 



ir(u>\y a [t'}\)^ P e [x-y a [t'}}, £>0 



Here, the prefactor 1/|^| 3 ^ 2 <C 1 exhibits Vainshtein screening of high frequency Galileon radiation. 
Counting powers of oj tells us that an 0[v 3 ]-accurate answer receives contributions from the I = 0, 1, 2 
terms: 
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(The m a d 2 y a /di /2 came from t = 1; the second group involving l/-y/w|y a | from 



v 7/2 „3™ 



(183) 
|&| > 1- 



2; and the final group 



from £ = with included.) Comparing equations (180) and (183), we find, as we did in the previous 
section, that the slope of the Galileon power spectrum depends on the size of At low frequencies 
the radiation is very similar to that of the non-interacting case, while at high frequencies, the spectrum 
acquires an overall suppression factor of 1/|£„| 3 / 2 . The phenomenology of Galileon radiation appears to 
be richer than its minimally coupled massless counterpart due to the existence of the additional length 
scale r v in the problem. Also note the presence of fractional powers of oj that cannot be associated 
with time derivatives - this aspect of the Galileon radiation has no analog in its non-interacting cousin 
nor in gravitational waves propagating on flat spacetime. 

Gravitational Waves Let us also record the power spectrum of GW emission. Arguments based 
on conservation of energy and the validity of Newton's third law leads us to infer that the lowest order 



that, in the r, r' <C r v limit, nonlinearities begin to render the solution offered by Q( ata - t,c ) invalid at distances closer to 
m than ^JmjMr' . 
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answer arises from the third time derivative of the quadrupole moment of the system. Misner, Thorne 
and Wheeler [35] equations 36.1 and 36.2 tell us 

4Hm) = J\ Q 'H • < 184) 

where is the Fourier transform of the triple time derivative of the quadrupole moment, 

WaM = J (vaAAVaAA ~ ~<%^]) ■ (185) 

As we have seen, the forces acting between compact bodies orbiting around a central mass M no longer 
obey Newton's third law if Galileons exist (%2 a m a d 2 y a /dt' 2 ^ 0) - this in turn means, eq. (184) may 



no longer be the leading order answer to the GW spectrum. We hope to return to this issue in the 
future. 

V. SUMMARY AND DISCUSSION 

Within the next decade or so, gravitational wave detectors are expected to begin hearing signals 
from astrophysical systems such as inspiralling compact binaries. It is therefore an appropriate time to 
explore the possibility that there could be emission of radiation due to additional degrees of freedom 
coming from the various modifications of General Relativity that has been proposed in the literature to 
date. Our work initiates such an investigation for Galileons, a class of scalar field theories that exhibit 
what is known as the Vainshtein mechanism, within the context where there is a large central mass M. 

We have constructed the Galileon retarded Green's function satisfying the linearized equations of 
motion about the background field of this central mass. The main results are described in section 



(III A). For the radiation problem, where the observer is situated at a very large distance from M and 



the source generating Galileon waves is located well within the Vainshtein radius of M, the primary 



results can be found in equations (44) through (48). We have also obtained the exact static Green's 



function, the WKB limit of the retarded Green's function, and the Green's function evaluated both 
deep inside and far outside the Vainshtein radius of M. 

We have used this radiative Green's function to obtain the frequency spectrum of Galileon radiation 



emitted from (acoustic) waves on the surface of the spherical mass M, described by equations (155), 



( p7l ), ( p8| ), and ( p9] ). The oower dissipated per unit solid angle per mode (£, m), in the non-relativistic 



limit, can be found in eq. ( 164 ) for the case where the r w -to-oscillation-time-scale ratio was much smaller 



than unity; and in (165) for that ratio very large. To illustrate the importance of the nonlinearities 
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of the Galileon interaction for the problem at hand, we have also calculated for comparison the power 
emitted if we replace the Galileon with a non-interacting massless scalar, with the relevant results given 



in equations ( 162 ) and ( 163 ) . 

A particularly interesting application of our results is to the radiation spectrum due to the motion 
of n point masses gravitationally bound to the central mass M . This is the dissipative aspect of the 
Galileon modified gravitational dynamics whose conservative portion we investigate in a separate paper 
[28]. We have focused on the non-relativistic limit, and have found the energy in Galileon radiation 



lost to infinity in two regimes. For small r„-to-orbital-time-scale ratios, this is given by eq. (178), and 



the lowest order answer is (180). For large r\,-to-orbital-time-scale ratios, the answer is eq. (181), with 



the lowest order result in eq. (183). For comparison, the non-interacting massless scalar result is eq. 



(177) and that of the quadrupole radiation formula for gravitational waves is eq. (184). 

In these radiative processes, we find that the non- interacting massless scalar and Galileon radiation 
are comparable in the non-relativistic, low frequency, and low multipole regime. Moreover, in this limit 
Galileon radiation is actually amplified relative to its non-interacting counterpart for higher multipoles 
- these findings are a direct consequence of the structure of the Galileon radial Green's function, 



which we have already highlighted in the discussion surrounding equations (49) and (50). In the 
high frequency limit, we confirm the anticipated Vainshtein screening of the Galileo radiation at low 



multipole orders; demonstrating it to be of O 



•3/2 



relative to its non-interacting counterpart. At 



high enough multipoles, however, high frequency Galileon radiation becomes enhanced relative to its 
non-interacting counterpart. Moreover, for the astrophysical n body system, where Newton's third 
law is obeyed between each compact body and M but not between the compact bodies themselves, the 
leading 0[v%\ terms in the radiation formulas, equations (180) and (183), are small (scaling as 0[m/M\) 
but non-zero, in contrast with the non-self-interacting scalar case, where ^ a m a d 2 y a /di /2 = 0. 

Having developed some quantitative understanding of the production of Galileon radiation in this 
paper, let us remark that, if Galileon waves exist, they are in principle detectable by GW detectors. 
About flat Minkowski spacetime, the sum of the graviton-matter and Galileon-matter coupling (eq. 
@) is 



M 



pi 



"fJlV 



h {eS) = h 



21177™ . 



(186) 



This implies that, if Galileons are present, ordinary matter experiences an effective weakly curved 
metric of the form 



h 



(cff) 



M, 



Pi 



(187) 
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and the tidal forces experienced by the arms of the interferometers of GW detectors would now be due 
to both the transverse-traceless graviton and the Galileon waves (p. 

Some of the calculations we have carried out have overlap with other recent work on the same topic 
[29]. In this paper, we have sought to understand the Galileon radiation spectrum in both the high 
and low frequency limits, and have found that the slope of the power spectrum of Galileon radiation 
should be a non-trivial function of £„, the ratio of the Vainshtein radius r v to the typical wavelength of 



the emitted waves. (For instance, the radiative limit of our retarded Green's function, namely, eq. (44) 
and the coefficients cl rad ^s described in (45) and (48), have very different ^-dependence for \£ v \ -C 1 



compared to 3> 1.) In [29], the focus was on the power loss from binary pulsar systems such as 
PSR B1913+16, and thus the authors carried out an analysis valid in the high frequency limit [38]. 12 
In the limit |£ t ,| 3> 1 we find that the Galileon power scales as |£u| -3//2 , in agreement with the results 
of [29]. 

We note that, a priori, caution is required in interpreting the results both in this paper and in [29] as 
describing a comparable mass binary pulsar system, because it is unclear if such a binary can be treated 
within the perturbative framework, given that the nonlinearities of Galileons are very important when 
the motion is taking place deep within the collective Vainshtein radius of the system itself. 13 However, 
the authors of [25] have explicitly verified the validity of their perturbative scheme, and will present it 
in an upcoming publication [381 139] . 

Future Work It would be worthwhile to convert the frequency space calculation in this paper 
to a real time one, in order to better understanding the physical meaning behind the fractional powers 
of angular frequency found in the Galileon emission spectrum. To this end, the identification of the 



correct contour prescriptions in the Fourier integral of eq. ( 39 ) would be necessary. We also have left 



12 As explained in [29] . for Galileons to be relevant for cosmology, it is often assumed that A ~ M p \Hq ~ l/(10 3 km), where 
Ho ~ 10 -33 eV is the current Hubble parameter. For binary pulars like PSR B1913+16, with masses on the order of a 
few solar masses, the corresponding Vainshtein radius is r v ~ C[10 3 ] light years. Because the period of typical binaries 
(~ 1/w) are of the order of a few hours, therefore |£„| ^> 1. 

In more detail, in [2§] the binary system, with masses Mi and M2, is modeled by adding and subtracting to the binaries' 
stress energies a monopole term with stress energy given by To[x] = (Mi + M2)S^[x]. The field generated by To is 
then used as a background, and the stress-energies of the pair of point masses themselves minus the stress energy of 
the central monopole (see their eq. (2.6)) are treated as perturbations. However, since, as in this paper, the linearized 
equation of motion about the background of the monopole are solved, the subtracted monopole really plays no role, and 
this scheme is really equivalent to the setup where there is one central mass M = Mi + M2 and two masses in orbit 
around it, one of mass Mi and the other M2. But since Mi and M2 are comparable in magnitude to Mi + M2, and there 
is no small dimensionless ratio one may use as an expansion parameter, it is not evident for the general binary problem 
that the Mi,2 are mere perturbations on top of the M. More quantitatively, recall the mode functions evaluated deep 
within the Vainshtein radius can be expressed in a separated form, reflecting the spherical symmetry of the background. 
In a comparable mass relativistic binary system, this spherical symmetry is completely absent. Thus, the separation of 
variables technique may not be useful in solving the general binary problem. 
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unexplored a large range of \£ v \, as we have only examined the extreme limits <C 1 and S> 1. 
Furthermore, since the nonlinearities of the Galileon theory play such a crucial role in its dynamics, we 
hope in the near future to go beyond the linear analysis about the background II due to M . We could 
compute, say within the Born approximation, the first correction arising from the cubic self-interactions 
m eq. ([6j) to the wave solutions we have obtained here, so as to better understand the domain of validity 
of the results in this paper. 

There are a number of interesting further directions for future work. In 4 dimensional flat spacetime, 
one should introduce the quartic and quintic Galileon terms and carry out an analogous analysis to 
that performed here. One may also wish to develop an understanding of the backreaction of the power 
loss on the motion of the n point masses. In this context, there are other processes one could consider 
- for example, one could carry out a calculation analogous to the one found in [36] and |37| for the 
gravitational case, in which one small mass scatters off M, producing Galileon bremsstrahlung radiation 
in the process. 
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Appendix A: Minimally Coupled Massless Scalar Radiation: Spacetime Calculation 



The central goal of this section is an alternate derivation of eq. (177), but (for technical simplicity) 
without the inclusion of gravitational interactions. We wish to compute the non-relativistic power 
emitted in non-interacting massless scalar radiation by the motion of n compact bodies. We will do so 
by finding the scalar field and its first derivatives generated by these point masses in flat spacetime; and 
from these gradients construct r 2 T 0r . However, we will first derive a expression in curved spacetime 
and proceed to specialize to Minkowski spacetime. The reason for doing so is that the Galileon ip 



field generated by the n body system we examined in section (IVB) is related to the problem of a 



non- inter acting scalar in the geometry given in eq. (128) 



Let the n masses be {m a \a = 1,2, ... ,n} and their spacetime locations {ya}- The massless scalar 
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theory in question is the curved spacetime generalization of the H-theory in eq. (144), namely 
with proper times 



(Al) 



a=l 



ds a = dt\ 



dt dt 



From the Hadamard form of the scalar Green's function in eq. (132), 

8[t-f 



with 



4tt 



U, 



{S[a]U x , x , + e[a]V XiX >) , 



(A2) 



(A3) 



(A4) 



(for an explanation of these symbols, see |34j ) . the solution to S is given by the integral 

n „ 

H ^ = J dSa0 [* - ^M] ( S Wx, Va Px,y a + ®Wx,y a ]V x ,y a ) , 

For a fixed x, cr Xt y a = constant defines the geodesic joining x (the observation point) to y a (the source 
point), provided this geodesic is unique (an assumption used when deriving the Hadamard form in eq. 
( |132 )). Thus it must be possible to invert <J x , ya for s a and vice versa. In particular, in the following, 
we will need 



der. 



ds a ds a 



V a/ a[x,x' = y a ] . 



and hence 



ds n 



da 



da 



ds a 



The null cone piece of H involves <5[cr X)3/a ] which we may then write 



8[s a - To] 



(A5) 



(A6) 



(A7) 



\da Xt y a /ds a \ 

where r a , the retarded time, is defined to be the proper time of the ath mass when it lies on the 
backward null cone of the observer at x, 



cr[x,y a [T a ]] = 0. 



(A8) 



The "backward" part of the requirement is reinforced by the step function [t — y^] in eq. ( A5 ) 
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The curved spacetime scalar field produced by the n point masses is therefore 



E 

a=l 



47rM p i I I \da Xt y a /ds a \ 



+ / ds'V[x,y a [s'} 



We may obtain the gradients of E by differentiating the integral representation in eq. (A5). 



V a E[x] = £ 4^ J dsa (v a a6'[a]U + 5[a]V a U + 5[a]V a aV + @[a]V a v\ . 



The b~'[o~\ term may be re- written, holding x fixed, 

ds a d S[s a - T a ] 



da X:Va ds a \da[s a = r a ]/ds a | ' 
which allows us to integrate by parts to obtain 



(A9) 



(A10) 



(All) 



V„C4, 



VaE[x] g 1,7.1/,,; [\ da x Jds a 



(da x ,yjds a ) 2 



P' 

a 

ds a 



+ Ti Pi ^ I -TyVA'ff^ + V K 'V\>(Tx,y a | + 



Id^^/dsal 3 I ds^ 



ds a ds a 



Idcr^/dsal 



(A12) 
+ / ds a V a V Xt y 



with U Xj3 j = y/ A XiX i. The primed derivatives are with respect to y a (the ath point mass location) and 
the unprimed ones with respect to x (the observer location). 

Minkowski spacetime In Minkowski spacetime, U = 1 and VJ7 = V = VV = 0, while the world 
function reads 



1 



^ T.T 1 



rj^ix - x'Y{x - x'Y 



(A13) 



so that its derivatives are 



V^ct^ = {x- y a )fj,, V^a x ,y a = (y a - x) M , V^V^^,^ = -77^, 



(A14) 



and 



der. 



x,y a 



(X ~ Va)i 



(A15) 



ds a ds a 

For convenience we shall use a dot to represent a derivative with respect to s a . The gradient of 
in Minkowski spacetime generated by n point masses is given by 



VH[a:] = - £ 



■ A' 

?7« 



^ 4vrM p i | (y£(y a - x) K ) 2 \y%(y a - x) K 



VaiVa- X)x + y a 



(A16) 
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where the retarded proper time r a is the solution of the equation 



t-ya[ T a} = \x-y a [r a ]\- 



(A17) 



Comparing equations (A12) and (A18), we see that in flat spacetime, the field detected by an observer 



(or felt by some test mass) at a; is a function of the positions of the n point masses evaluated at 
the appropriate retarded times - that is how long it takes the signal to reach the observer from the 
source. Whereas in curved spacetime, where Huygens' principle no longer holds, the problem of motion 



is a richer one, because the proper time integral involving the tail term V x , Va in eq. (A12) receives 
contributions from the point masses' entire past histories. (The electromagnetic counterpart to this 



statement can be found in |40| . in which the and F^ v analogues of equations (A9) and (A12) are 
derived.) 

We are now ready to determine the power radiated by the n point masses moving in flat spacetime. 
Let us assume the motion of these n bodies is confined within some finite spatial volume containing 
the center of the spatial coordinate system; this will certainly be true when these n bodies are bound 
together by their mutual gravity, which is the central theme of this paper. We wish to extract the piece 
of H that represents radiation. As already explained earlier, in Minkowski spacetime, the radiative 



piece of H is the 1/r piece. From the result in eq. (A18), we see that this can be identified by counting 
powers of (x — y a ); this isolates the acceleration term. If we let the observer lie at some very large 



radius, we deduce y^iVa — x)) 



ValVa + Va- {x - y a ) -» r(-y% + y a ■ x) and y%(y a - x) x 



\x - ValVa + Va-{x- y a ) -> r(-y° a + y a -x). Recalling the notation dy%/ds a = y% and d 2 y%/dsl = y r a , 



oo 



(l,x) ^ m a VaVafS {±,xf 



47rr ^M p i |y0_y> a .£|3 



(A18) 



and the power radiated to infinity per unit solid angle r 2 T 0r 

dE 



2G N / >A y*r] a p(l,xf 



r 2 V r HV'H is 

2 



i/i 3 



(A19) 



To take the non-relativistic limit, we first recall that d/ds a = (1 - (dy a /di') 2 ) _1/2 d/dt / , where the 



retarded time t' satisfies eq. (A17), i.e. t' = t 



y a [t']\. That retarded time depends on the 



trajectory y a means, upon carrying out the retarded time derivatives with respect to t' in eq. (A19), 



we still need to Taylor expand every time dependent expression in powers of |y a [t']|/r, because the 
latter expansion introduces further time derivatives. It is at this point that terms containing the same 
number of time derivatives (now evaluated at the approximate retarded time t' = t — r) are considered 
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to be of the same order in the non-relativistic expansion. Up to 0[v^] we have 

) ■ (A20) 



t'=t-r/ 



Appendix B: Green's Function of Linear Second Order Ordinary Differential Equation 

In this section we review the construction of the Green's function for linear second order ordi- 
nary differential equations (ODEs) in terms of their homogeneous solutions. Consider the differential 
operator 

d 2 d 

^z = P2[z]-^+Pi[z]—+ Po [z], (Bl) 

where po,i,2[-z] are smooth functions of z. The homogeneous solution / to the corresponding linear 
second order ODE satisfies D z f[z] = 0, and the associated symmetric Green's function obeys 

® x G[z, z'\ = ® z ,G[z, z'\ = X[z]5[z - z'\ = X[z']5[z - z'\ . (B2) 

To solve this, first assume z > z f , so that S[z — z'\ = 0. The problem is thus reduced to solving 
the homogeneous equation D z G[z,z'] = "£) z iG[z,z'] = 0. Let /i^I-z] be a pair of linearly independent 
homogeneous solutions, i.e. 

®zfi[z] = Q z f 2 [z]=0, Wr w [/i,/ 2 ]^0, (B3) 
where the Wronskian Wr is defined to be 

Wr (a) [/i,/ 2 ] = fi[z}(f2)'[z} - {h)'[z\h[z\. (B4) 
Because D z G[z, z'\ = 0, we must have 

G[z> z']=a I f 1 [z], 1 = 1,2 (B5) 
where the o^s are ^-independent. Similarly, D z >G[z,z'] = means 

^ = ^[4 J = 1,2, (B6) 
where A> is now a 2 x 2 matrix of z, z'-independent constants. The same argument would hand us, for 



z < z', 



G[z < z'\ = A l lh[z\h[Ai I, J = 1,2 (B7) 
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If G[z, z'\ were not continuous at z = z' its second derivative with respect to z or z' there, and hence 



the second line of eq. (B2 ), would involve a derivative of S[z— z'\. That implies we may assume G[z, z'] 



is continuous at z = z'. A}?fi[z]fj[z] = A l > fj[z]fj[z] imposes the conditions 



All _ All a22 _ a22 



21 



A% + A > 



21 



A< + A< 



(B8) 
(B9) 



We now integrate D z G[z, z'\ = X[z]5[z — z'\ around the neighborhood of z = z', applying integrating-by- 
parts as many times necessary to shift all the derivatives acting on G[z, z'\ onto the P2,l- By continuity, 
thepo term, \piG] z z_ z ',_ , [(j? 2 )' G] z z. z ',_ and the remaining integral involving G (with no derivatives acting 
it) drop out - assuming pi 2 are smooth - leaving us with 

=*'+ 



XW]=P2[Z'} 



dG[z, z' 
dz 



(BIO) 



Employing the continuity conditions in equations (B8 ), one would find the and Arr do not contribute 



to X[z']. We may use (B9) to eliminate, say Al^, and find (BIO) becomes 



X[z'} = - P2 [z'](A 21 - A^)\Nr {zl) [h, f 2 ] 



(Bll) 



Because any "rotation" of the pair /i^, i.e. the pair {qi = fj\l = 1, 2} for any 2x2 invertible Q, is 
an equally valid pair of linearly independent homogeneous solutions, without loss of generality we may 



choose A^ = A^} — 1, such that now 



X[z} = - P2 [z]\Nr {z) [f 1 ,f 2 



(B12) 



Therefore, the general solution to eq. (B2) is 



G[z,z')=Cf 1 [z > ]f 2 [z < ] - (1 -C)/i[*<]/ 2 [*>] 
+ C 11 f 1 [z]f 1 [z'} + C 22 f 2 [z}f2[z'], 



(B13) 



where z> (or z<) is the greater (or smaller) of the pair (z,z ! ), and C, C 11 and C 22 are arbitrary 
constants. These coefficients will be fixed by the boundary conditions of the given physical problem. 
By using the differential equation obeyed by the fi 2 , one may readily show that 



d Wr (z) [/ 1 ,/ 2 ] = -^M\A/r {2) [/ 1 ,/ 2 ] 



(B14) 



dz'" (Z)lJ1,J ^ 1 p 2 [z] 

This in turn means the Wronskian of two linearly independent solutions, and hence the measure X[z], 
can be solved up to an overall constant, without first solving for the homogeneous solutions. Recalling 
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eq. (B12), we gather 



AN = X P2p]exp 



dz 



//PIP 



p 2 [z" 



(B15) 



where x is t ne constant. In solving for G[z, z'], we will thus first choose a value for x, and use this 



choice to fix the normalization of the product of the solutions /i[^>]/2[^<] in eq. (B13). 

We wish to emphasize, because we are constructing a symmetric Green's function, one that is a 
Green's function with respect to both variables z and z' , we have just seen that one does not have a 
choice in picking the measure X[z], but rather X[z] is fixed up to an overall numerical constant by p 2 [z] 



in eq. (Bl) and the Wronskian of any two linearly independent homogeneous solutions. 

To summarize, once a pair of homogeneous solutions f\ 2 are known, the symmetric Green's function 



has the general solution given in equation (B13). The measure X[z] multiplying the ^-functions in eq. 



(B2) is given by eq. (B15), and the overall constant x there needs to be chosen. In a given problem - 
for this paper it is the solution of ge[€, £'] i n e Q- Q - the analog of C, C n , and C 22 (or the A 1<2 and 
-81,2) will be fixed by appropriate regularity and boundary conditions. Because X[z] has been computed, 
the overall normalization of the products /i[£>]/2[#<] is determined by the Wronskian condition in eq. 



(B12). 



5 -Function Measure Let us conclude this section by justifying the (rr') 1 measure on the right 



hand side of eq. (31 ). Equation (B15 ) informs us that we can determine this measure up to a constant, 



by integrating the ratio of — 2e^/^ to — e 2 . This may be achieved by using the explicit expressions in 

f< 2e 3 



p8| through (l30|) yielding 



exp 



£"e2[£",&] 

This immediately tells us that the measure multiplying the (^-functions in eq. (31) is 



Xe2[£,6u]exp 



x_ 
£ 2 ' 



(Bi6) 



(B17) 



where x is a constant. Far away from the central mass M, using equations (19) and (20) to keep only 



the most dominant terms in 61,2,3 (equations (28) through (30)), the left hand side of eq. (31) yields, 
as expected, the flat spacetime minimally coupled massless scalar wave equation 



ri^d^d u G[x,x'} = —^[t - t']d[r - r']5[cos6 - cos8'}5[$ - (f>'} . 



(B18) 



Therefore choosing x = 1 amounts to adhering to the usual convention. 
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